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1.  INTRODUCTION  AND  EXECUTIVE  SUMMARY 


This  document  is  the  final  report  (Item  No.  0001  AC)  for  SDIO  contract  number  N000I4- 
90-C-0145,  "Versatile  High  Performance  Holographic  Optical  Coatings."  The  contract  has  been 
monitored  by  ONR.  The  purpose  of  the  work  was  to  investigate  the  potential  for  fabricating 
holographic  coatings  (HCs)  with  nearly  arbitrary  reflectivity  profiles.  The  desired  reflectivity 
spectrum  of  the  HC  could,  for  example,  be  the  optimum  matched  filter  for  detecting  a  specified 
target  against  a  bright  high  clutter  background.  Once  the  desired  spectral  profile  of  the  filter  was 
calculated  using  well  known  matched  filter  theory,  holographic  techniques  would  be  used  to 
produce  the  filter.  The  analysis  presented  in  this  report  has  shown  that  sucn  a  procedure  will 
work  and  allows  the  filters  to  be  fabricated  much  more  easily  than  would  be  the  case  for 
conventional  multi-layer  dielectric  coatings,  deposited  one  layer  at  a  time. 

Section  2.  contains  the  background  material  which  examines  how  filters  can  be  used  in 
target  detection.  Simple  examples  include  the  very  narrow-band  filters  employed  in  laser  radar 
and  somewhat  wider  band  IR  filters  used  to  detect  missile  plumes  against  an  Earth  background. 
In  laser  radar,  the  filters  block  all  radiation  wavelengths  from  reaching  the  detector  except  the 
narrow-band  return  laser  wavelengths.  Missile  plume  detection  in  the  IR,  against  a  high  clutter 
Earth  background,  is  best  accomplished  at  the  strong  atmospheric  water  absorption  line  centered 
at  2.8  pm  or  the  strong  C07  absorption  line  centered  at  4.3  pm.  The  Earth  background  is 
several  orders  of  magnitude  dimmer  at  these  wavelengths  (filter  widths  ~0.2  pm),  allowing 
missile  detection  once  the  target  breaks  cloud  cover.  The  filters  for  these  applications  can  easily 
be  made  with  conventional  multi-layer  coatings,  but  HCs  offer  an  advantage  even  for  these 
simple  cases.  These  advantages  include: 

•  Simplified  fabrication,  especially  for  higher  performance  filters  where  a  fine 
spectral  structure  is  desired  within  the  narrow-band. 


Higher  damage  thresholds,  both  from  nuclear  events  and  from  hostile  laser 
irradiation. 


In  most  of  the  HC  applications  investigated  in  this  contract,  the  HC  is  used  in  a  reflector 
mode.  The  target  and  background  radiation  passing  through  the  camera  beam  train  is  reflected 
from  the  HC  and  onto  the  detector.  The  spectral  filtering  is  thus  performed  through  reflection 
and  not  transmission  of  the  HC.  The  Fresnel  reflection  of  the  HC  can  easily  be  eliminated  by 
recording  the  HC  at  an  angle  so  that  the  interference  planes  of  the  coating  are  tilted  slightly  with 
respect  to  the  surface  of  the  HC. 

The  real  advantage  of  the  HC  occurs  for  more  complex  applications  than  narrow-band 
filters.  It  is  very  difficult,  for  example,  to  detect  missile  plumes  in  the  visible  against  the  Earth 
background,  which  is  especially  intense  in  the  visible.  In  a  matched  filter,  the  filter  will  allow 
wavelengths  where  the  target  is  especially  strong,  and  the  background  less  intense,  to  be 
transmitted  to  the  focal  plane.  It  will  also  attenuate  those  wavelengths  where  the  target  is  dimmer 
and  the  background  is  brighter  from  reaching  the  focal  plane.  Several  examples  are  outlined  in 
section  2  of  this  report  which  indicate  that  a  matched  filter  could  be  used  to  detect  missile  plumes 
in  the  visible,  contrary  to  the  generally  held  belief  of  many  "IR  detector  people."  There  are 
significant  advantages  to  detecting  missiles  in  the  visible  or  near  IR.  These  include: 

•  Cryogenic  cooling  is  not  required  for  silicon  detectors  used  in  the  visible,  thus 
greatly  reducing  the  cost,  complexity,  and  weight  of  such  a  system,  compared  to 
IR  focal  plane  arrays. 

•  The  silicon  detectors  represent  a  much  more  mature  technology  than  HdCdTe  or 
InSb  detectors.  They  are  also  much  less  expensive,  partly  because  of  economics 
of  scale.  Silicon  CCDs  are  mass  produced  for  a  commercial  market. 

Section  2  concludes  with  examples  of  how  holographically  produced  filters  could  be  used 
commercially.  The  use  of  these  filters  could  help  identify  specific  minerals  and  vegetation  from 
a  surveillance  satellite. 

In  section  3,  two  basic  concepts  are  examined  for  fabricating  the  general  purpose  HCs. 
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The  analysis  for  both  concepts  is  performed  using  the  coupled  mode  theory.  The  final  equations 
which  determine  the  coating’s  reflectivity  spectrum  in  terms  of  exposure  wavelength,  angle,  and 
exposure  time,  must  be  solved  numerically.  Before  this  is  done,  we  review  a  number  of 
candidate  hologram  recording  materials  in  section  4  of  this  report.  The  most  attractive  material, 
because  of  its  high  diffraction  efficiency  (large  index  of  refraction  modulation),  high  resolution 
(>10,000  lines/mm),  low  scattering,  and  moderate  exposure  levels  is  dichromated  gelatin 
(DCG).  Section  4  reviews  the  properties  of  DCG  and  in  particular  how  an  unusual  property  of 
DCG  can  be  exploited  to  simplify  the  filter  fabrication  for  the  "Concept  2"  approach,  described 
in  section  3. 

Section  5  contains  the  results  of  our  numerical  analysis.  We  considered  "Concept  2"  HCs 
designed  for  two  military  applications: 

•  Detection  of  missile  plumes  in  the  UV  (3000  to  4000A). 

•  Detection  of  missile  plumes  in  the  visible  against  an  Earth  background. 

The  section  5  analysis  indicates  that  HCs  can  be  made  easily  that  accurately  represent  the 
ideal  spectral  profiles  needed  for  the  UV  and  visible  plume  detection.  Some  of  the  key  results 
of  this  program  have  been  submitted  for  publication  and  presentation  at  the  1991  SPIE 
International  Symposium  on  Optical  Applied  Science  and  Engineering  (21-26  July  1991,  San 
Diego,  CA;  Conference  No.  15).  A  draft  copy  of  the  paper  is  given  in  Appendix  5B. 
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2.  BACKGROUND 


Before  discussing  the  novel  approaches  for  fabricating  optical  filters  which  have  been 
examined  during  Phase  I  of  this  program,  we  will  review  some  of  the  potential  uses  of  these 
filters  in  the  fields  of  surveillance,  target  detection,  and  tracking.  The  primary  emphasis  during 
Phases  1  and  2  of  this  program  will  be  in  the  development  of  the  holographic  filters  for  military 
applications,  especially  for  SDI  applications.  These  applications  will  be  discussed  in  section  2.1. 
During  Phase  3  of  this  work,  as  is  often  the  case  in  SBIR  programs,  a  greater  emphasis  will  be 
placed  on  commercial  uses  of  the  holographic  filters.  These  uses  are  examined  in  section  2.2. 

2.1.  MILITARY  APPLICATIONS  FOR  MATCHED  HOLOGRAPHIC  FILTERS 

2.1.1.  MISSILE  LAUNCH  DETECTION  IN  THE  IR 

One  of  the  most  well  known  military  uses  of  optical  filters  is  in  missile  launch  detection 
in  the  IR.  The  bright  plume  of  a  missile  in  the  IR  can  be  used  to  detect  a  missile  from  space, 
even  with  a  bright  Earth  background,  when  viewed  with  a  moderately  narrow  band  optical  filter 
(AA  -  0.2  pm)  centered  at  one  of  the  absorption  lines  of  the  atmosphere.  For  example,  the 
Earth  background  intensity  drops  several  orders  of  magnitude  around  2.8  pm  and  4.3  pm,  due 
to  strong  water  and  carbon  dioxide  absorption  lines,  respectively.  An  optical  filter  centered  at 
one  of  these  lines  will  therefore  block  most  of  the  Earth  background  radiance,  but  transmit  the 
radiation  from  a  plume,  once  the  missile  breaks  cloud  level  and  is  above  most  of  the  attenuating 
atmosphere. 

Of  course,  these  optical  filters  can  and  have  been  made  with  conventional  multi-layer 
deposition  techniques  (e.g.  physical  vapor  deposition),  but  holographic  filters  offer  a  number  of 
advantages.  These  include: 

•  Potentially  narrower  transmission  lines  and  simplified  fabrication, 
especially  for  advanced  filters  whose  detailed  spectral  properties  could  be 


2-1 


matched  to  the  atmospheric  transmission  in  order  to  achieve  even  higher 
signal-to-clutter  ratios  (SCR). 

•  Higher  damage  thresholds,  both  from  nuclear  events  and  from  hostile  laser 
irradiation. 

The  principal  investigator  for  this  SBIR  Phase  I  program,  Dr.  Jeffrey  B.  Shellan, 
experimentally  demonstrated  the  high  laser  damage  thresholds  of  holographic  coatings  under  a 
subcontract  to  SAIC  (No.  DAAM01-89-R-0065),  which  was  funded  by  SDIO  in  1989. 

2.1.2.  LASER  RADAR 

Another  military  (as  well  as  non-military)  application  of  optical  filters  is  in  laser  radar 
systems.  The  laser  radar  can  be  used  for  target  ranging,  estimating  target  velocities  (from  the 
Doppler  shifted  return),  and  for  background  clutter  suppression.  For  example,  once  a  missile’s 
bright  plume  is  detected,  a  future  space  interceptor,  such  as  Brilliant  Pebbles  may  use  laser  radar 
to  intercept  the  missile.  The  ladar  could  initially  be  used  to  obtain  ranging  information  needed 
for  calculating  the  optimum  intercept  trajectory.  As  the  two  vehicles  converge,  a  "hard  body 
handoff"  is  required— the  interceptor  must  lock  onto  the  missile  body  itself,  rather  than  the  much 
brighter  plume.  One  approach  for  ensuring  that  this  handoff  is  made  is  to  use  laser  radar  for 
actually  imaging  the  hard  body.  A  relatively  large  region  near  the  tip  of  the  plume  is  illuminated 
by  the  laser  and  the  scattered  laser  light  is  detected  by  a  CCD  covered  with  a  very  narrow  band 
filter.  The  filter,  which  is  centered  on  the  laser  line,  blocks  nearly  all  the  radiation  from  the 
plume  and  Earth  background  and  detects  the  scattered  laser  light.  For  a  liquid  fueled  missile, 
very  little  laser  radiation  is  scattered  from  the  plume  itself,  so  centroid  tracking  can  be  used  for 
hard  body  intercept.  For  solid  fueled  missiles,  however,  the  fuel  particles  in  the  plume  may 
actually  scatter  as  much  laser  radiation  as  the  missile.  In  this  case,  the  target  must  be  imaged 
using  the  focal  plane  array  and  an  algorithm  used  which  can  identify  the  hard  body  silhouette. 

Of  course,  the  optical  filter  needed  for  the  above  mission  could  be  fabricated  by 
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conventional  means,  but  holographic  filters  offer  some  advantages.  In  addition  to  those  listed 
at  the  end  of  section  2.1.1,  narrow  line  HCs  can  be  made  with  larger  fields-of-view  (FOV)  than 
standard  coatings.  This  important  point  will  be  discussed  in  the  next  section. 

2.1.3.  FRAUNHOFER  LINE  SURVEILLANCE  -  SINGLE  FOCAL  PLANE 


To  first  order,  the  solar  spectrum  can  be  approximated  by  a  5900°K  black  body.  This 
spectrum  does  have  a  number  of  very  narrow  (~  1  -  10A)  spectral  holes,  however,  caused  by 
solar  atmospheric  absorption.  These  regions  of  attenuated  solar  intensity  are  referred  to  as  the 
Fraunhofer  lines  and  the  author  of  this  report  has  investigated  exploiting  this  effect  in  target 
detection.  This  work  was  done  as  part  of  a  two  year  JBST  contract  from  Lawrence  Livermore 
National  Laboratory. 

Some  of  the  most  well  known  Fraunhofer  lines  are  the  H-a  line  (6562. 8A),  the  Mg  II 
lines  (2795. 4A  and  2802. 3A),  and  the  Ca  II  lines  (8498. lA,  8542. lA,  and  8662. 2A).  At  the 
center  of  the  above  six  lines,  the  solar  intensity  is  attenuated  by  84%,  90%,  90%,  70%,  81%, 
and  79%,  respectively.  "Single  focal  plane  Fraunhofer  line  surveillance”  is  realized  by  covering 
a  focal  plane  with  a  narrow  band  filter  centered  at  one  of  the  Fraunhofer  lines.  Any  background 
clutter  produced  by  reflected  sunlight  will  thus  be  attenuated  by  70  to  90%.  The  self  radiating 
targets,  such  as  missile  plumes,  however,  are  not  effected,  since  their  radiation  does  not  originate 
from  the  sun.  This  filter  will  thus  increase  the  SCR  by  a  factor  of  5  to  10.(1)  Larger  apertures 
are  required  for  this  surveillance  scheme  to  keep  the  SNR  at  an  acceptable  level,  since  this 
approach  improves  the  SCR  at  the  expense  of  the  SNR. 

It  is  doubtful  that  the  target  detection  capability  (SCR  and  SNR)  of  a  "Fraunhofer 
detector"  could  ever  equal  that  of  the  conventional  IR  approach  described  in  section  2.1.1.  There 


(,)SCR  increases  of  two  to  three  orders  of  magnitude  can  be  realized  using  "dual  focal 
plane  Fraunhofer  line  surveillance,"  to  be  described  in  section  2.1.4. 
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is,  however,  one  very  important  advantage  of  the  Fraunhofer  concept  over  the  IR  approach  -  the 
detection  can  be  done  with  visible  wavelength  detectors.  The  potential  advantages  of  using 
visible  wavelength  Si  detectors  is  so  great  that  LLNL  has  spent  a  considerable  amount  of  time 
investigating  a  number  of  these  novel  schemes  for  their  Brilliant  Pebbles  program.  The 
important  advantages  of  visible  wavelength  detectors  over  IR  detectors  include: 

•  Lower  cost  (economics  of  scale  due  to  mass  production  of  visible  wavelength 
CCDs  for  commercial  use) 

•  Lighter  weight 

•  Well  developed  technology  for  large  focal  plane  arrays  with  high  quantum 
efficiency 

•  Cryogenic  cooling  is  not  required 

Unfortunately,  the  Earth  background  is  much  brighter  in  the  visible  and  typical  targets 
are  much  dimmer  in  the  visible,  than  in  the  IR.  Thus  novel  techniques  such  as  "Fraunhofer 
detection"  will  be  needed  if  this  economical  target  detection  approach  is  to  be  realized. 

A  key  difficulty  which  JBS  Technology,  Inc.  encountered  in  designing  the  "Fraunhofer 
detector,"  as  well  as  some  other  classified  approaches  to  visible  wavelength  target  detection,  was 
the  incompatibility  of  narrow  line  filters  and  wide  fields-of-view.  For  laser  radar  filters,  narrow 
receiver  FOVs  are  often  acceptable  since  the  outgoing  laser  beam  must  be  quite  narrow  in  its 
angular  extent  to  ensure  an  adequate  return  signal.  For  passive  warning  and  detection  systems, 
however,  wide  FOVs  are  often  needed. 

Figure  2.1  illustrates  the  spectral  transmittance  of  a  narrow  line  interference  filter, 
designed  for  transmission  at  0.500  pm  when  illuminated  at  normal  incidence.  As  the  filter 
illumination  angle  is  increased,  the  peak  transmittance  shifts  to  shorter  wavelengths  and  broadens. 
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This  characteristic  of  interference  filters  would  severely  restrict  a  system’s  FOV,  which  used  a 
narrow  band  filter.  Figure  2.2  is  a  plot  of  a  filter’s  FOV  as  a  function  of  its  percentage  half 
width.  The  half  width  of  the  H-a  Fraunhofer  line  is  about  2k  or  about  0.03%  of  the  wavelength 
(2/6562).  Figure  2.2  indicates  that  the  full  cone  FOV  of  this  system  would  only  be  about  4°  for 
a  filter  with  an  effective  index  of  3. 

One  method  for  increasing  the  surveillance  camera’s  FOV  would  be  to  put  the  narrow 
band  filter  on  a  curved  surface  as  illustrated  in  figure  2.3.  This  approach  was  investigated  by 
the  author  of  this  report  several  years  r  jo.  hut  we  were  unable  to  find  a  single  vendor  able  to 
deposit  a  multi-layer  ("conventional")  coating,  to  the  required  layer  thickness  control,  over  a 
large  curved  substrate.  This  coating  could  easily  be  made,  however,  using  holographic  methods. 
A  photosensitive  material,  such  as  dichromated  gelatin  (DCG),  would  be  deposited  on  the  curved 
glass  substrate.  A  converging  and  diverging  beam  could  then  be  interfered  in  the  DCG  to 
produce  the  desired  interference  filter. 

2.1.4  FRAUNHOFER  LINE  SURVEILLANCE  -  DUAL  FOCAL  PLANE 


In  section  2. 1.3  we  saw  that  the  background  clutter  could  be  reduced  by  70%  to  90%  by 
the  use  of  narrow  band  filters  centered  at  Fraunhofer  lines.  It  has  been  shown  by  JBS 
Technologies,  Inc.  that  the  use  of  a  pair  of  focal  planes  covered  by  twu  different  narrow  band 
filters,  one  centered  at  a  Fraunhofer  line  and  the  other  filter  shifted  several  angstroms  off  of  the 
Fraunhofer  line,  can  increase  the  SCR  by  at  least  three  orders  of  magnitude!  A  similar 
improvement  in  SCR  could  be  realized  by  using  atmospheric  absorption  lines  in  the  near  IR.  such 
as  the  water  absorption  line  at  0.92  pm  (visible  wavelength  Si  CCDs  can  be  used  up  to  ~  1 
pm). 


Figure  2.4  is  a  schematic  diagram  of  the  hardware  needed  for  this  concept.  The  signals 
from  the  two  co-registered  focal  planes,  which  measure  very  closely  spaced  spectral  bands,  can 
be  combined  in  an  optimum  manner  which  nearly  cancels  the  background  clutter  completely. 
Some  details  involved  in  this  signal  processing  are  discussed  in  Appendix  2A.  Appendix  2A  also 
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presents  an  example  of  the  power  of  this  scheme  using  real  Earth  scene  data,  as  previously 
prepared  by  JBST  under  a  LLNL  contract. 

It  is  commonly  believed  that  boosters  cannot  be  detected  at  visible  wavelengths  from 
space  with  a  bright  Earth  background.  Although,  for  reasonably  sized  camera  apertures,  the 
booster  plume  radiant  intensity  is  high  enough  to  generate  an  acceptable  SNR,U)  the  background 
clutter  is  too  large  to  allow  detection  by  conventional  means-the  SCR  is  too  low. 

Figure  2.5  is  a  plot  of  the  detection  range  as  a  function  of  camera  lens  diameter  for  a 
third  stage  booster/1 2’  as  seen  at  0.656  pm  (the  Fraunhofer  H-a  line).  We  assume  that  both  a 
SNR  and  SCR  of  5  are  needed  for  detection.  The  plot  indicates  that  if  no  filters  are  used,  the 
detection  range  is  only  80  km,  so  the  booster  could  not  be  detected  at  cloud  break,  as  desired, 
from  a  surveillance  satellite  at  400  km.  Note  that  the  detection  range  of  80  km  is  independent 
of  the  camera  lens  diameter,  since  it  is  SCR  limited  and  not  SNR  limited.  Larger  lens  diameters 
allow  more  target  photons  to  be  collected,  but  more  background  photons  are  also  collected  so  the 
detection  range  cannot  be  increased.  The  use  of  a  narrow  band  filter  centered  at  the  H-a 
Fraunhofer  line  and  a  single  focal  plane  can  increase  the  detection  range  to  about  120  km,  since 
the  background  intensity  can  be  preferentially  reduced  with  this  filter.  This  detection  range  is 
still  too  limited  to  be  of  interest. 

If  a  dual  focal  plane  with  complementary  narrow  band  holographic  filters  is  used, 
however,  the  detection  range  can  be  extended  to  several  hundred  kilometers.  A  10  cm  aperture 
provides  a  large  enough  signal  to  allow  detection  at  550  km.  Since  nearly  all  the  background 
clutter  can  be  suppressed  with  this  scheme,  the  camera's  performance  is  limited  by  SNR  and  not 


(1) That  is,  readout  noise,  detector  noise,  shot  noise,  etc.  are  ail  low  enough  relative  to  the 
target  signal. 

(2) The  booster  type  will  not  be  given  in  order  to  keep  this  report  unclassified. 
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the  SCR.  This  is  an  important  conclusion  since  visible  wavelength  detection  of  boosters  offers 
some  advantages  over  IR  detection,  as  outlined  in  section  2.1.3. 

2.1.5.  WIPE-BAND  MATCHED  FILTER  FOR  BOOSTER  DETECTION 


Most  of  the  SDI  related  filter  designs  discussed  thus  far  have  been  narrow  band  filters. 
For  the  detection  of  dim  targets  against  a  fairly  uniform  background,  it  may  be  desirable  to  use 
a  wider  band  filter  so  that  more  target  photons  can  be  transmitted  to  the  focal  plane.  A  "matched 
filter”  has  a  spectral  profile  which  blocks  transmission  of  radiation  to  the  focal  plane  at 
wavelengths  where  the  background  is  bright  and  the  target  is  dim.  It  has  a  high  optical 
throughput  at  wavelengths  where  the  target  is  brighter  and  the  background  is  dim.  If  the  spectral 
profile  of  the  background  and  target  are  known  apriori,  then  a  filter  with  the  'deal  spectral  profile 
can  be  designed.  As  a  final  example  of  a  military  application  for  holographic  filters,  we  will 
consider  the  design  of  such  a  broad  band  filter  for  the  detection  of  a  solid  fueled  missile  against 
an  Earth  background. 

A  solid  fueled  missile  plume  can  be  modeled,  to  first  order,  as  a  2350°K  black  body. 
The  Earth  background,  in  the  visible,  will  be  modeled  as  a  5900°K  black  body  (reflected  solar 
spectrum)  which  has  been  modified  by  the  spectral  transmission  of  the  Earth’s  atmosphere. 

Figure  2.6  is  a  normalized  plot  of  the  black  body  curves  for  a  5900°K  radiator,  whose 
spectral  peak  is  at  0.5  pm,  and  a  2350°K  radiator,  with  a  peak  at  ~  1.2  pm.  Figure  2.7  is  a 
plot  of  the  Earth  background  radiance,  found  by  taking  the  product  of  the  spectral  profile  of  the 
5900° K  black  body  and  the  atmospheric  transmission  (double  pass),  as  calculated  using  the 
Lowtran  7  atmospheric  model.  Two  different  curves  are  plotted  in  figure  2.7  -  one  for  the 
Lowtran  Mid-Latitude  Summer  atmospheric  model  and  one  for  the  Sub-Arctic  Winter 
atmospheric  model. 

If  the  Earth  background  profile,  as  shown  in  figure  2.7,  is  combined  with  the  2350°K 
plume  profile,  as  plotted  in  figure  2.6,  an  optimum  filter  spectral  profile  can  be  calculated.  This 
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matched  filter  is  the  ideal  filter  for  detecting  the  taiget.  The  reflectivity  profile  of  the  desired 
holographic  filter  is  shown  in  figure  2.8.  We  have  designed  the  filter  to  operate  in  a  reflectivity 
mode,  so  the  detector  receives  radiation  reflected  by  the  hologram  (rather  than  transmitted).  Two 
different  reflectivity  profiles  are  indicated  in  figure  2.8,  corresponding  to  the  two  different 
atmospheric  models.  The  figure  shows  that  it  may  be  desirable  to  use  a  filter  wheel  containing 
several  holographic  filters.  The  best  filter,  depending  on  atmospheric  conditions,  could  then  be 
rotated  in  place  in  front  of  the  detector. 

A  more  attractive  application  of  a  holographic  matched  filter  for  plume  detection  can  be 
realized  when  the  targe*  has  some  detailed  spectral  structure.  Examples  of  this,  in  the  UV 
regime,  are  illustrated  in  figure  2.9.  A  holographic  match*''*  filter  cold  be  fabricated  which 
would  pass  a  large  portion  of  the  spectrum  where  the  missue  plume  is  intense  and  block  the 
wavelengths  between  the  intensity  spikes  shown  in  figure  2.9. 


2.2  COMMERCIAL  APPLICATIONS  FOR  MATCHED  HOLOGRAPHIC  FILTERS 

Most  of  he  Phase  II  efforts  will  be  directed  toward  developing  HCs  which  have  a  direct 
military  application.  This  technology  can  also  be  used  in  commercial  applications,  such  as 
improved  identification  of  specific  vegetation  or  minerals  from  satellite  cameras.  In  fact  airborne 
imaging  spectrometers  have  already  been  used  to  identify  various  types  of  vegetation  (and  even 
the  health  of  the  vegetation)  and  numerous  minerals--all  from  the  "spectral  fingerprint"  of  the 
specific  material  being  searched  for.  This  technology  could  eventually  be  used  for  worldwide 
surveys  to  predict  crop  yields,  for  monitoring  illegally  grown  crops  (coca,  marijuana,  the  opium 
poppy),  and  to  conduct  rapid  geological  surveys  for  mining  and  soil  type  delineation.  It  could 
be  used  to  detect  vegetation  stress,  to  monitor  tropical  forest  changes,  to  perform  forest  structure 
analysis  and  forest  biomass  estimates,  and  even  to  measure  plankton  type  and  concentration  (to 
predict  the  location  of  fish).  This  could  all  be  done  from  orbiting  satellites. 

The  AIS  (airborne  imaging  spectrometers)  hardware  being  developed  by  JPL  for  these 
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measurements  uses  a  dispersing  grating  to  direct  various  wavelength  bands  to  the  N  rows  of  a 
2-D  focal  plane  array.  Each  row  of  the  CCD  thus  images  the  same  1-D  strip  of  the  scene  being 
viewed,  to  obtain  a  spectrally  resolved  1-D  view  of  the  target  area.  The  2-D  spatial  image  is 
generated  by  "push  broom  imaging"  of  the  Earth  scene  passing  below  the  camera.  Since  most 
materials  have  a  characteristic  spectral  fingerprint,  this  technology  represents  a  major  advance 
in  remote  sensing.  The  concept  is  illustrated  in  figure  2. 10.  Figure  2. 1 1  and  2.12  illustrate  this 
spectral  fingerprint  for  several  different  plants  and  minerals,  as  measured  with  a  spectrometer. 

In  the  imaging  spectrometer  (IS)  concept,  a  separate  image  is  generated  at  numerous 
wavelengths  (over  200  separate  wavelengths  in  the  advanced  JPL  devices)  and  this  wealth  of  data 
is  then  processed  to  identify  material  in  the  scene.  Commercial  applications  of  holographic 
matched  filter  technology,  which  we  hope  to  develop  in  Phase  III  of  this  program,  will 
complement  the  IS  systems.  Rather  than  generating  hundreds  of  images  of  the  same  scene,  as 
viewed  at  different  wavelengths,  only  one,  or  at  most  a  handful,  of  images  are  taken.  Rather 
than  using  very  narrow  band  filters  or,  equivalently,  a  dispersive  grating,  filters  with  complex 
spectral  properties  could  be  used  in  the  holographic  matched  filter  concept.  A  holographic  filter 
would  be  designed  to  highlight  a  specific  material  in  the  scene.  The  approach  is  not  as  powerful 
as  IS,  but  is  much  simpler,  requires  far  less  computer  power,  and  may  be  more  useful  for  rapid 
surveys  of  large  areas  where  it  is  desirable  to  use  the  entire  2-D  focal  plane  for  imaging,  and  not 
just  one  row. 

A  filter  wheel  placed  in  front  of  the  CCD  could  contain  special  holographic  matched 
filters  designed  to  highlight  areas  of  the  scene  with  a  desired  spectral  characteristic.  For 
example,  the  DEA  may  wish  to  mount  a  coca,  marijuana,  and  opium  holographic  filter  on  their 
camera,  while  the  Agriculture  Department  might  want  space  reserved  for  a  com  or  soybean 
filter!  These  holographic  filters  have  the  property  that  they  allow  radiation  to  reach  the  focal 
plane  at  wavelengths  for  which  the  material  of  interest  radiates  strongest.  The  filters  block  the 
transmission  at  wavelengths  for  which  the  material  does  not  radiate  strongly.  As  was  the  case 
for  the  military  targets  discussed  in  section  2.1,  much  higher  SCR  enhancement  will  be  realized 
if  a  dual  focal  plane  configuration  is  used  with  complementary  holographic  matched  filters. 
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We  believe  that  many  of  the  filters  described  in  the  last  two  sections,  especially  ones  with 
complex  spectral  properties,  can  best  be  fabricated  holographically.  Now  that  the  motivation  for 
the  development  of  these  complex  filters  has  been  presented,  we  will  describe  some  approaches 
for  fabricating  them  and  analyze  their  performance. 
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Figure  2. 1  Angular  properties  of  all-dielectric  interference  filters.  Measured  variation  with  angle 
of  incidence  of  the  spectral  transmittance  of  a  typical  commercial  interference  filter. 
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Figure  2.3  Fabricating  a  narrow  band  filter  over  a  curved  surface  can  be  used  to  increase 
FOV  of  a  surveillance  camera. 
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Figure  2.4  Schematic  diagram  of  two  focal  plane  spectral  signal  processing  concept 
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THE  COMBINATION  OF  AN  H-ALPHA  FRAUNHOFER  FILTER  AND  SIGNAL  PROCESSING 
ALLOWS  THE  DETECTION  OF  A  THIRD  STAGE  BOOSTER  AT  VISIBLE  WAVELENGTHS 
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Normalized  black  body  distribution  F(  X)  (meters-1 ) 


Wavelength  in  nm 


Figure  2.6  Plot  of  the  Normalized  Black  Body  Curves  for  2350°K  and  5900°K 
Radiators 
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OPTIMUM  MATCHED  FILTER  REFLECTIVITY  AS  A  FUNCTION  OF 
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Figure  2.9  Examples  of  UV  Spectra  From  Storable  Propellant  Rocket  Plumes. 
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IN  100-200  SPECTRAL  BANDS, 
INHERENTLY  REGISTERED 


WAVELENGTH,  fj. m 
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Figure  2.11  Fluorescence  Spectra  of  Plants  of  Four  Major  Plant  Types  Provide 
Information  Regarding  Type  of  Plant  Stress  and  Identification  of  Specific 
Plant  Types.111  (Active  IS) 
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^Alexa:  F.H.  Goetz,  et  al,  "Imaging  Spectrometry  for  Earth  Remote  Sensing,"  Science  22S,  1147  (1985 

Figure  2.12  Examples  of  Short  IR  Spectra  of  Twelve  Minerals111  (Passive  IS) 
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APPENDIX  2A.  MATCHED  FILTER  SIGNAL  PROCESSING 


2A.1  THEORY 


This  Appendix  will  summarize  some  of  the  mathematical  details  which  support  our 
contentions  in  section  2. 1 .4  that  SCR  enhancements  of  several  orders  of  magnitude  are  achievable 
if  complementary  optical  filters  cover  co-registered  focal  plane  arrays.  This  spectral  sampling 
and  signal  processing  scheme  does  require  that  the  spectral  bands  be  judiciously  chosen. 
Examples  include  the  Fraunhofer  lines  (section  2.1.4)  or  Earth  atmospheric  absorption  lines. 

We  will  define  the  signal  received  in  pixel  number  i,  focal  plane  number  j,  as  en.  The 
signal  e;j  depends  not  only  on  the  scene  being  observed,  but  also  on  the  optical  transmission 
property  of  the  filter  covering  focal  plane  j.  The  energy  measured.  eij7  consists  of  a  possible 
target  signal  t(j  and  a  background  signal  b-: 


e.. 

u 


=  b,i 


+  t. 


We  wish  to  form  a  linear  sum  of  the  signals  from  pixel  number  i  from  each  focal  plane 
and  will  define  this  processed  signal  as  s;: 


(t) 

s  =  s  +  s 


n(b> 
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where: 


0 


j-i 


3. 


s<b)  = 


j=i 


a.  b.. 

j  'j 


4. 


The  sums  in  equations  3  and  4  are  over  the  two  focal  planes  in  our  system.  The 
coefficients  aj  to  be  used  in  the  signal  processor  are  chosen  so  that  the  variance  of  the 
background  signal  over  the  field  of  pixels  is  as  small  as  possible,  while  at  the  same  time,  the 
processed  target  signal  Sj(1)  is  as  large  as  possible.  This  is  a  classical  matched  filter  problem,  and 
the  details  involved  in  solving  the  problem  will  not  be  presented  here.  They  were  presented  by 
Dr.  Shellan  in  reference  1.  The  result  is  given  by: 


B*'t  5. 

a  =  — — 

~  tT  B"‘t 


In  equation  5,  the  filter  coefficients  a^  are  expressed  by  the  column  vector  a  whose 
components  are  the  desired  coefficients  a,.  The  column  vector  t  represents  the  total  target  signal 
measured  in  each  of  the  two  focal  planes.  The  superscripts  "T"  and  "-1"  denote  transpose  and 
matrix  inverse,  respectively.  The  matrix  B  is  found  from  background  spectral  correlation 
properties,  and  the  j,  k  component  of  this  matrix,  Bjk,  is  related  to  the  bliin  coefficients  by: 

1  p 

B  •  1  Y  b  b ,  - 

jk.  ij  ik 

P  i  =  l 


1  y  b.. 
p  h  ,J 


i  * 

iy  b 
p  h  lk 


6. 


The  integer  p  in  the  summation  in  equation  6  represents  the  number  of  pixels  in  our  field 
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of  interest.  The  output  target  signal  from  our  signal  processor  E  s[°  ,  divided  by  the  rms 

..  > 

background  signal,  is  defined  as  the  signal-to-clutter  ratio  of  the  signal  processor.  It  can  be 
shown(1>  to  be  given  by: 

7. 

SCR  =  (tTB-,t)1/2 

Now  that  we  have  outlined  some  of  the  mathematics  involved  in  this  problem,  we  will 
explain  in  more  physical  terms  how  these  results  will  be  used  to  design  a  system  best  able  to 
detect  targets.  The  parameters  that  determine  the  signal-to-clutter  ratio  (SCR  in  equation  7)  are 
the  total  target  energy  measured  by  each  focal  plane  and  the  Earth  background  spectral  statistics 
as  given  by  the  matrix  elements  Bjk. 

The  matrix  elements  Bjk  simply  represent  various  spectral/spatial  correlation  properties 
of  the  background  scene.  Analysis  indicates  that  for  the  novel  sampling/signal  processing  scheme 
proposed,  the  output  SCR  is  proportional  to  at/AA,  where  the  variables  a,  t,  and  A  A  represent 
the  atmospheric  absorption,  target  signal,  and  separation  of  the  two  sampled  wavelengths, 
respectively.  Thus,  we  expect  that  the  SCR  can  be  driven  to  very  high  levels  by  exploiting 
narrow  atmospheric  absorption  lines  which  in  turn  allow  the  use  of  closely  spaced  sampling,  A  A. 
Fraunhofer  lines  may  allow  even  smaller  A  A  values  to  be  used. 

2A.2  CONCEPT  DEMONSTRATION  AND  NUMERICAL  ANALYSIS 

Figures  2A.la  and  2A.lb  show  a  pair  of  coregistered  and  radiance-calibrated  images 

•‘^The  variable  o  actually  represents  the  differences  in  the  atmospheric  absorption  between 
the  target  and  more  distant  background.  If  the  sampling  is  done  with  a  Fraunhofer  line,  then  a 
represents  the  strength  of  the  Fraunhofer  line  dip. 
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collected  over  the  Australian  desert  by  the  TIMS  sensor  in  two  adjacent  LWIR  wavebands 
separated  by  A  k  =  400  nm  (center-to-center).  The  measured  spectral  covariance  matrix  for  these 
scenes  (denoted  as  g  in  section  2. A.  1)  is 


B11  B12 

47.16  49.50 

B  = 

~ 

B12  B22 

49.50  52.38 

This  reveals  that  the  clutter  has  nearly  the  same  average  power  in  each  band  and  an  extremely 

high  correlation  coefficient  of  p12  =  Bp  /^Bn  B,2  =  0.996  between  bands.  If  a  target  response 

were  present  with  roughly  equal  magnitude  t  in  bands  1  and  2,  very  little  could  be  gained  by 
combining  the  bands  on  a  pixel-by-pixel  basis.  For  example,  with  t  =  7,  the  SCR  in  either  band 
is  on  the  order  of  unity,  the  signal  vector  is  tT  =  [7  7],  and  from  equation  7,  the  largest  SCR 
available  from  combining  both  bands  is 


SCR  =  ( tT  B‘ 1 1 ) 1/2  -  0.96 

which  is  far  too  low  for  reliable  detection 


9. 


However,  if  the  clutter  response  in  one  of  the  two  bands  (say,  Band  1  in  figure  2A.la) 
is  synthetically  attenuated  to  simulate  the  effect  of  selective  background  absorption,  the  situation 
changes  significantly.  Figures  2A.2a  and  2A.2b  show  grey-scale  plots  for  the  central  64  x  64 
portion  of  the  two  images  in  figure  2A.  1 .  If  the  apparent  pixel  radiances  in  the  Band  1  subimage 
of  figure  2A.2a  are  scaled  by  the  factor  o  =  0.5  but  those  in  Band  2  are  left  unchanged,  then 
the  measured  spectral  correlation  matrix  becomes 


Bll  B12 

11.79 

24.75 

B  = 

' 

BI2  B22 

24.75 

52.38 
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Repeating  the  SCR  calculation  using  equation  10  with  the  same  signal  tT  =  [7  7],  we  now  find 
that 


SNR  =  (tTB-'t)1/2  -  12 


11. 


To  demonstrate  that  such  a  signal  is  now  readily  detectable,  the  unresolved  additive  target 
response  shown  in  figure  2A.2c  is  synthetically  injected  into  both  the  attenuated  Band  1  subimage 
and  the  Band  2  subimage  with  a  magnitude  of  7  units.  Figure  2A.2d  shows  the  results  of  taking 
the  optimum  linear  combination  of  the  two  resulting  scenes,  as  defined  by  equation  5  (the  output 
image  is  plotted  on  the  same  grey-scale  as  figures  2A.2a  and  2A.2b  to  facilitate  a  direct 
comparison). 

In  analysis  conducted  by  the  principal  investigator,  expressions  for  the  SCR  gain  of  the 
proposed  concept  were  developed  as  a  function  of  a,  A k,  and  a  set  of  physical  parameters  which 
are  related  to  the  measured  spectral  correlation  of  the  background  in  the  two  bands.0'  For 
example,  using  the  data  in  the  two  TIMS  images  in  figure  2A.I  to  estimate  the  mode! 
parameters,  the  achievable  SCR  gain  against  this  terrain  background  can  be  locally  parameterized 
by  a  and  A  A  as  shown  in  figure  2A.3.  This  family  of  curves  clearly  demonstrates  the  benefits 
of  selective  absorption  and  small  wavelength  differentials  for  a  specific  case.  Note  the  relatively 
high  SCR  achievable  even  with  absorption  coefficients  as  small  as  5%  for  small  A  k  values. 


(1)Jeffrey  B.  Shellan,  "Clutter  Suppression  Through  Differences  in  Target  and  Background 
Spectral  Content,"  J.B.S.  Technologies,  Inc.,  report  to  Lawrence  Livermore  National 
Laboratory,  Contract  B056075  (October  12,  1988). 
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Figure  2A. la  TIMS  Band  1  LWIR  Image  of  the  Figure  2A.  1b  TIMS  Band  2  LWIR  Image  of  the 

Australian  Desert  (10  m  ground  resolution)  Australian  Desert  (10  m  ground  resolution) 


I 

I 


Figure  2A.2c  Point  target  response 
Injected  Into  the  Band  1  and  Band  2 
subimages 


Figure  2A.2d  Results  of  optimum  linear  processing  of 
Bands  1  and  2  with  Injected  target.  (Band  1  background 
clutter  attenuated  by  a  =  0.5  to  simulate  selective 
absorption.) 
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RELATIVE  SCR  GAIN 


Figure  2A.3  SCR  gain  vs.  wavelength  differential  (AA.)  for  various  absorption 
coefficients  a.  Gain  is  measured  relative  to  SCR  obtained  by  iinear 
processing  of  two  bands  outside  absorption  region. 
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3.  EXAMINATION  OF  APPROACHES  FOR  FABRICATING  HOLOGRAPHIC  MATCHED 
FILTERS 


3.1  NARROW-BAND  FILTERS 


The  simplest  HC  to  fabricate  is  the  narrow-band  filter  which  reflects  radiation  over  a 
narrow  spectral  band  and  transmits  nearly  all  the  radiation  outside  of  this  band.  This  kind  of 
filter  has  been  fabricated  for  years  and  is  reviewed  here  simply  as  a  means  to  introduce  the 
complex  filters  to  be  developed  in  this  program. 

The  HC  is  fabricated  by  interfering  two  counter  propagating  plane  waves  in  a  photo¬ 
sensitive  material  such  as  dichromated  gelatin  (DCG).  This  is  usually  accomplished  by  placing 
a  reflector  plate  behind  the  hologram  being  recorded  so  that  the  incident  wave  is  reflected,  thus 
producing  the  desired  standing  wave  pattern.  This  geometry  is  illustrated  in  figure  3.1.  The 
intensity  modulation  will,  n  turn,  produce  an  index  of  refraction  modulation  in  the  developed 
hologram.  The  perod  -  i  the  index  modulation,  Ag,  is  equal  to  Ar/(2n0),  for  a  hologram  recorded 
with  normally  incident  light  (A.  =  wavelength  of  recording  radiation  in  vacuum,  nu  =  average 
index  of  refraction  of  photosensitive  layer).  This  index  of  refraction  period  is  the  period  which 
will  produce  the  greatest  reflectivity  when  the  hologram  is  subsequently  re-illuminated  with 
radiation  at  wavelength  kr  (normal  incidence). 

For  hologram  recording  at  oblique  incident  angles,  as  shown  in  figure  3.2,  the  fringe 
spacing,  Ag,  will  be  greater  than  for  normal  incident  recordings.  For  a  reflection  hologram 
recorded  at  an  incident  angle  4>r,  the  period  of  the  index  modulation  will  be: 


2  (n2a  -  sin2  4>r) 


1/2 


3.1 


The  above  grating  will  only  reflect  strongly  at  wavelength  Ar,  under  subsequent  re¬ 


s' 1 


illumination,  if  the  "playback"  incident  angle  is  also  4>r.  For  normal  incidence  operation,  the 
grating  period  given  in  equation  3.1  will  reflect  most  strongly  when  illuminated  by  radiation  of 
wavelength  Ap  given  by: 


k  = 

p 


n  k 

o  r 


( n;  -  sin2  4>r ) 


1/2 
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It  is  thus  straightforward  to  record  a  hologram  at  wavelength  kT  which  can  strongly  reflect 
at  a  longer  wavelength,  k  ,  when  the  hologram  is  later  illuminated  at  normal  incidence.  If  the 
DCG  is  sequentially  exposed  at  a  number  of  incident  angles,  <J>,,  <t>2,  ...  4>n,  by  radiation  at  a 
single  wavelength,  A.r,  then,  after  development,  the  resulting  hologram  will  strongly  reflect  at  a 
number  of  discrete  wavelengths,  k{,  k2,  ...  A.n,  determined  from  equation  3.2.  This  effect  has 
been  demonstrated  by  many  groups  and,  in  fact,  is  the  basis  for  some  filters  designed  to  protect 
space  surveillance  sensors  from  blinding  hostile  laser  radiation.  It  has  also  been  used  to  fabricate 
laser  goggles  which  can  provide  the  user  with  eye  protection  from  several  common  laboratory 
lasers.  We  will  now  show  how  this  concept  can  be  used  to  fabricate  holographic  filters  with 
arbitrary  complex  reflectivity  profiles,  such  as  those  needed  for  the  "matched  filters"  described 
in  section  2. 

3.2  CONCEPT  1:  CONTINUOUS  EXPOSURE  WITH  PHASE  CONTROL 

3.2.1  DESCRIPTION  OF  CONCEPT  1  HOLOGRAM  FABRICATION 


In  section  3. 1  it  was  pointed  out  that  if  the  holographic  recording  medium  is  sequentially 
exposed  to  recording  beams  at  a  number  of  different  incident  angles,  the  index  of  refraction 
modulation  of  the  developed  hologram  will  consist  of  a  superposition  of  gratings,  each  with  a 
different  period.  This  structure  will  then  reflect  strongly  at  a  number  of  different  wavelengths. 
For  normal  incidence,  the  hologram  will  reflect  strongly  at  wavelengths  equal  to  2n0  A?lp),  where 
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is  the  period  of  the  grating  within  the  recording  medium. 


The  basic  concept  to  be  investigated  in  this  program  is  whether  a  hologram  with  an 
arbitrary  reflectivity  profile  can  be  fabricated  by  sequentially  recording  gratings  with  the  desired 
period.  To  fabricate  a  hologram  with  a  continuous  reflectivity  spectral  profile,  R(A),  the 
hologram  could  be  recorded  by  exposing  the  DCG  (or  other  photosensitive  material)  to  a 
continuum  of  incident  angles,  before  the  DCG  is  actually  developed.  The  strength  of  the 
reflectivity  desired  at  a  specific  wavelength  would  be  controlled  by  modulating  the  intensity  of 
the  laser  used  to  record  the  hologram,  as  the  angle  of  incidence  of  the  recording  beam  was 
changed.  The  geometry  for  such  a  system  is  illustrated  in  figure  3.3.  The  above  approach  for 
fabricating  a  hologram  with  any  desired  reflectivity  profile  will  be  referred  to  as  "Concept  1." 

In  the  following  section  we  will  analyze  the  expected  performance  of  the  novel 
holographic  filter  fabricated  by  the  procedure  described  above.  Obviously  the  key  issue  to 
investigate  is  how  the  continuum  of  "off  resonance"  Fourier  grating  components  effects  the 
reflectivity  of  the  beam.  Although  it  has  been  experimentally  demonstrated  that  recording  a 
"reasonable"  number  of  different  grating  periods  in  one  hologram  can  produce  a  structure  that 
will  then  reflect  at  each  of  the  corresponding  radiation  wavelengths,  this  has  not  been 
demonstrated  in  the  limit  of  an  infinite  number  of  different  grating  periods— that  is,  for  the  case 
of  a  continuum  of  grating  periods. 

3.2.2.  CONCEPT  1  ANALYSIS 

3. 2.2.1  RELATIONSHIP  BETWEEN  REFLECTANCE  PROFILE.  r(k).  AND  INDEX 
PROFILE,  nlz) 

This  problem  will  be  analyzed  using  the  coupled  mode  theory,  as  described  in  the 
references  of  Yariv,(1)  Yeh,(2)  Shellan,(J)  and  Kogelnik.(4)  We  will  briefly  summarize  this  theory. 
The  propagation  of  radiation  is  described  by  the  wave  equation,  which  for  an  electric  field  with 
time  dependence  exp(-itot)  reduces  to: 
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+  k2  n2  E  =  0 


where: 


E(z)  -  electric  field  in  the  holographic  coating 

z  ■  coordinate  of  a  point  within  the  coating,  as  measured  normal  to  the 

coating/air  interface  (z  =  0  is  the  coordinate  at  the  coating  surface) 
k  ■  radiation  wave  number  in  vacuum  =  to/c,  where  w  is  the  radian 

frequency  of  the  radiation  and  c  is  the  speed  of  light 
n  =  n(z)  =  index  of  refraction  of  the  HC 

If  the  hologram  were  recorded  with  light  at  a  single  angle  of  incidence,  then  the  index 
of  refraction  could  be  represented  by  a  sinusoidal  modulation  of  the  form: 


i  ,  .  2  ~  2  7X 

n*(z)  =  nc  +  2non,cos  —  z  +  y 

\ 


where: 


(for  n,  <  <  n()) 
average  index 

amplitude  of  index  modulation 

grating  period 

phase  factor  for  grating 
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In  our  case,  however,  the  hologram  will  be  created  by  sequentially  exposing  the  DCG  to 
radiation  from  a  continuum  of  angles,  as  illustrated  by  the  hardware  in  figure  3.3.  In  this  case 
the  index  of  refraction  of  the  DCG  after  development  must  be  represented  by  a  continuum  of 
grating  periods: 


(max) 

g 

n2  (z)  ■  no  +  2no  J  n,  (Kg)  exp(iKgz)dKg 

j,(min) 


(max) 

Kg 

+  2no  [  n_  (Kg)  exp(-i  Kgz)dKg 

K(min) 


where: 


K 


g 


=  grating  vector 


Kg(m,n)  and  Kg(max)  are  the  minimum  and  maximum  values  of  the  grating  vector 
which  make  up  the  complex  hologram  profile. 

n,  and  n  are  the  Fourier  coefficients  which  describe  the  hologram  profile. 


The  third  expression  needed  to  derive  the  coupled  mode  equations  is  the  equation  for  the 
form  of  the  electric  field.  The  field  E(z)  will  be  written  as: 


E(z)  ■  F(z)  exp  ( i  no  kz )  +  B(z)  exp  (  -i  no  kz  ) 


3.5 
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In  equation  3.5,  F(z)  represents  the  amplitude  of  the  incident  wave  as  it  travels  within  the 
hologram  and  B(z)  represents  the  amplitude  of  the  reflected  wave  within  the  hologram.  In 
coupled  mode  analysis,  we  use  the  fact  that  B  and  F  are  relatively  slowly  varying  functions  of 
z,  so  that: 


d2F 

dz2 


<  <  k2  n^  F 
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d2B 

dz2 


<  <  k2n 


B 
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If  we  combine  equations  3.1,  3.4,  and  3.5  and  ignore  the  second  derivatives  of  B  and  F, 
we  Find: 


2  i  n  k  —  exp  (in  kz )  -2  in  k  —  exp  (-in  kz  ) 
°  dz  F  0  °  dz  F  0 


+  2  no  k2  [  F  exp  ( i  no  kz  )  +  B  exp  (  -i  nu  kz )  ]  • 


(max) 

g 

•  J  n.(K)  exp  ( i  Kg  z )  dKg 

(min) 

g 


2no  k2  [  F  exp  ( i  no  kz )  +  B  exp  (  -i  no  kz )  ]  •  3  15 


..(max) 

g 

J  n_  (K,)  exp  (  -i  Kgz)dKg  =  0 

K<min) 

g 
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Since  F  and  B  vary  relatively  slowly  with  z  and  the  terms  exp(i  nQ  k  z)  and  exp(-i  nQ  kz) 
vary  rapidly  with  z,  equation  3. 15  can  be  broken  up  into  two  different  equations.  This  is  the  key 
step  in  the  derivation  of  the  coupled  mode  equations  and  is  discussed  in  detail  in  any  one  of  the 
four  references  given  above.  The  result  is: 


rr  (max) 

^■g 

^  -ikB  J  n+(Kg)  exp  [  i(  Kg  -  2n0  k)z  ]  dKg  =  0 

Z  j^(min) 
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wlmax' 

Kg 

^  +ikF  J  n_(Kg)  exp  [  -  i(  Kg  -  2nQ  k)z  ]  dKg  =  0 
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Equations  3.16  and  3.17  cannot  be  solved  analytically  for  the  general  case,  but  we  can 
find  an  approximate  analytic  solution  in  the  low  reflectivity  limit.  This  case  is  referred  to  as  the 
undepleted  incident  wave  approximation,  since  we  replace  F(z)  with  its  value  at  z  =  0  in 
equation  3.17.  Using  F(z)  =  F(0)  =  1,  equation  3.17  can  be  rewritten  as: 


j^(max) 

9  3  18 

^  ~  -ik  J  n_  (Kg)  exp  [  -i(  Kg  -  2n0  k)z  ]  dKg 

j,  (min) 

N 

If  we  integrate  each  side  of  equation  3.18  over  z,  from  z  =  0toz  =  L(L  =  holographic 
coating  thickness),  and  use  the  boundary  conditions: 

B(L)  =  0 

B(0)  =  r  =  coating  reflectance. 

we  can  relate  the  coating’s  reflectance  profile  r(k)  to  the  Fourier  component  h_ .  The  result  is: 
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(max) 

Kg 

r(k)  -  i  k  L  J  fl.  (Kg)  exp  [i(2nok  -  Kg)  L/2  ] 

„(min) 

Kg 


•  sine  [  ( 2no  k  -  Kg )  L/2  ]  d  Kg 


{  .  sin 

sine  u  *  — 
(  u 


3.19 


For  large  values  of  L  (thick  holographic  coatings),  the  sine  term  in  the  integrand  of 
equation  3.19  acts  like  a  delta  function.  Thus  if: 


3.20 

(2nok  -  Kgmin) )  L/2  >  >  1 

3.21 

and  ( 2  no  k  -  Kg™’0 )  L  "l  <  <  1, 
equation  3.19  can  be  approximated  as: 


r(k)  -  i2kn_(2nok)  J  sine  u  ellJ 


du 


3.22 


n.  (2nok)  -  —  r(k) 

71  Iv 


3.23 


Equation  3.23  shows  that  the  Fourier  component  n  (Kg)  is  proportional  to  the  coating's 

reflectance  evaluated  at  the  wavenumber  k  =  Kg/(2n0).  If  we  combine  equations  3.4  and  3.23, 
and  use  the  fact  that  n2(z)  is  real,  we  can  find  what  the  coating’s  index  of  refraction  profile  must 
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be  in  terms  of  its  reflectance  spectrum.  The  result  is: 


7  /  \  2 

n-(z)  -  nc  + 


4no 


‘‘max 


i  J  _  r  (k)  exp  ( -i  2no  kz )dk  +  c  c 


3.24 


where: 


kmin 

kmax 

"cc 


Kg(mi">/(2n0) 

Kg(max)/(2n0) 

represents  complex  conjugate 


For  relatively  small  index  modulations,  equation  3.24  can  be  rewritten  as: 


2  n 

n(z)  -  no  +  - 

71 


i  j"  I  r(k)  exp( -i2nokz)dk  +  c  c 
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3.2.2  2  CALCULATION  OF  HOLOGRAM  EXPOSURE  LEVELS 


In  this  initial  analysis  we  will  assume  that  the  index  of  refraction  modulation  of  the 
developed  hologram  is  proportional  to  the  exposure  of  the  DCG  during  the  recording  process. 
That  is,  non-linear  terms,  which  can  be  especially  important  if  the  DCG  is  over  exposed,  will 
not  be  considered  here.  If  6n  represents  the  change  in  the  index  from  its  unexposed  value,  then 


Sn  =  c  H 

O 


3.26 
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where: 

H  ■  exposure  of  DCG  (Joules/cm2) 

cQ  =  proportionality  constant  between  H  and  5n 

If  ID  is  the  intensity  of  the  exposure  beam  used  to  create  the  hologram,  then  the  standing 
wave  pattern  formed  in  the  hologram  during  exposure  will  have  an  intensity  ranging  from  0  to 
4I0  (This  assumes  that  the  two  counterpropagating  beams  forming  the  standing  wave  pattern  each 
have  intensity  IQ).  If  H  represents  the  DCG  exposure  that  would  result  from  only  one  of  the  two 
counterpropagating  beams,  then  the  total  peak-to-valley  index  modulation  produced  by  the 
standing  electric  field  is  4c0H. 

For  fabrication  "Concept  1,"  illustrated  in  figure  3.3  and  described  in  section  3.2.1,  the 
DCG  is  exposed  through  a  range  of  angles  by  rotating  the  mirror  indicated  in  figure  3.3.  The 
rotating  mirror  will  turn  at  a  constant  rate  during  the  exposure  and  the  laser  output  will  be 
modulated  so  that  the  exposure  level  is  a  function  of  the  angle  of  incidence  of  the  laser.  This 
will  control  the  amplitude  of  the  various  grating  periods,  which  in  turn  determines  the  reflectivity 
spectral  profile  of  the  hologram  after  development.  We  will  now  calculate  what  the  intensity  of 
the  exposing  laser  beam  should  be  as  a  function  of  incident  angle,  in  terms  of  c0  and  the  desired 
hologram  reflectivity  profile. 

If  H(0)d0  represents  the  differential  hologram  exposure  due  to  a  single  beam  as  it  swings 
through  the  angular  range  0  to  0  +  d0,  then  the  modulation  of  the  square  of  the  index  of 
refraction.  An2,  can  be  expressed  as: 


where: 

n„ 

Y(0) 


3.27 

An2  -  4  no  co  J  H  (0)  cos  ^  Kg  (0)  z  +  y  (0)  I  d0 
=  average  DCG  index  of  refraction 

=  a  phase  term  which  represents  the  relative  phasing  of  the  Fourier 
components  which  form  the  index  of  refraction  profile 
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The  integration  in  equation  3.27  is  over  the  range  of  recording  angles.  The  term  Kg(0) 
is  equal  to  2u/Ag(0),  where  Ag(0)  is  the  period  of  the  grating  formed  by  a  recording  beam 
incident  at  angle  0.  If  the  recording  wavelength  is  (in  vacuum)  and  a  superstrate  of  index  ns 
covers  the  DCG  during  the  exposure,  then  Kg  is  given  by: 


K  (0)  =  *1  [nl-n;  sin2  0  ] 1/2 


3.28 


If  we  compare  equations  3.24  and  3.27,  we  can  find  an  expression  for  H(0)  in  terms  of 
r(k).  The  result  is: 


4noco  J  H(0)  cos  [  Kg  (0)z  +  Y(6)]d0 

4  n„  ri 

—  i  -  r(k)  exp  (  -i  2n  kz )  dk  +  c  c 

7t  J  k 


3.29 


If  we  express  cos(Kg  z  +  y)  as  -  { exp  [  i  (  Kg  z  +  y  )  ]  +  c  c  }  and  equate  integrands 


on  each  side  of  equation  3.29,  we  find: 


2  n  l  i  fir- 

H  (0)  =  -  exp(iy)  i  r(k)  i 

Ti  c.  k  d0 
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The  HC  reflectivity.  R(k),  is  related  to  the  reflectance.  r(k),  by: 


3.31 

r(k)  =  v'R(k)  exp  [-i(J(k)] 


where  P  =  reflectance  phase  shift 


Since  H(0)  must  be  real  and  positive,  we  must  have: 


i  exp  [i(y  -  p  )]  -  -1 
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and  thus: 


H  (0) 


-2  n 

o 

nek 

O 


^R(k)  — 
v  d0 
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The  wave  vector  k  =  k(0)  in  equation  3.33  is  equal  to  Kg(0)/(2no),  where  K„(0)  is  given 
by  equation  3.28.  We  can  use  these  relationships  to  calculate  dk/d0.  Also,  if  the  "rotating 
mirror"  in  figure  3.3  turns  at  an  angular  velocity  of  Q  (radians/ second),  the  exposure  variable 
H(0)  is  related  to  the  incident  intensity  at  angle  0,  1(0),  by 


1(0)  =  QH(0) 


3.34 


Combining  equations  3.28,  3.33,  3.34  and  k  =  Kg/(2n0)  yields  our  desired  result: 


1(0) 


2  n  n.2  Q 

O  s 

"Co 


sin  0  cos  0 
(  n;  -  ns2  sin2  0  ) 


1  1/2 


R 


2  n 
n  A. 


(  nj  -  n:  sin2  0  ) 


1/2 


3.35 


Equation  3.35  indicates  exactly  what  the  recording  laser  intensity  should  be,  at  each 
exposure  angle  0,  in  order  to  create  a  hologram  with  the  desired  reflectivity  spectrum  R(k).  The 
required  intensity  is  proportional  to  Q/c(,  and  increases  as  R  increases,  as  expected. 


3-12 


3. 2.2.3  "INDEX  STACKING" 


Equation  3.35  indicates  that  a  coating  with  an  arbitrary  reflectivity  spectrum  can  be 
created  by  properly  modulating  the  recording  beam  as  the  hologram  exposure  angle  is  changed. 
Equation  3.25  gives  the  required  index  profile  as  a  function  of  the  desired  reflectance.  We  will 
now  apply  equation  3.25  to  a  specific  case  to  see  if  things  are  as  simple  as  they  appear! 

We  will  consider,  as  an  example,  the  case  of  a  coating  with  a  uniform  reflectivity  of  10% 
over  the  visible  wavelength  spectrum  (0.4  to  0.7  pm)  and  zero  reflectivity  elsewhere.  If  we  take  p  =  ~ 

so  that  r(k)  =  i  \/R(k)  -  +0.316  i  (equation  3.31)  and  evaluate  the  index  of  refraction  n(z)  at 
the  hologram  surface  (z  =  0),  then  equation  3.25  predicts: 

n(0)  =  no 

-  n(0)  =  no 

n(0)  =  no[ 1  +  0.225  ] 


1-1  |  r(k)  |  fn  (  kinix/kmin ) 

71 

1  +  1  X  0.316  fn  ( 0.7/0. 4  ) 


For  DCG,  the  average  background  index,  n0,  is  about  1.5,  so  equation  3.36  indicates  that 
the  index  at  the  surface  of  the  hologram  must  be  equal  to  1.84,  or  0.34  higher  than  the  average 
DCG  index.  This  level  of  index  modulation,  however,  is  impossible  to  achieve  in  DCG  or  any 
other  holographic  recording  material!  Experiments  have  shown  that  the  highest  peak-to-vallev 
index  modulation  achievable  in  DCG  is  about  0.2.  The  linear  relationship  between  DCG 
exposure  and  index  modulation  shown  in  equation  3.26  does  not,  of  course,  hold  for  arbitrarily 
large  exposures.  A  more  accurate  model  for  the  relationship  between  6n  and  H  is  a  non-linear 
one: 
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6n  =  6n 


1  -  exp 


'  -c.H' 


6n 


3.37 


For  small  exposure  levels,  equation  3.37  reduces  to  the  linear  result  given  in  equation 
3.26,  but  as  H  increases  the  index  change  saturates  at  the  value  &nmax. 


Achieving  a  10%  reflectivity  at  a  "single  wavelength"  requires  only  a  small  index 
modulation  for  a  holographic  coating  with  a  reasonable  thickness.  For  example,  the  peak-to- 
valley  sinusoidal  index  modulation  would  only  need  to  be  about  5  x  10'3  to  achieve  a  reflectivity 
of  10%  at  A.  =  0.5  pm  for  a  20  pm  thick  coating.  It  was  the  broad-band  reflectivity  requirement 
(0.4  to  0.7  pm)  which  drove  the  index  modulation  (at  the  coordinate  z  =  0)  up  to  0.34  for  the 
example  given  above.  In  order  to  produce  the  broad-band  spectrum,  the  index  of  refraction 
profile  must  contain  Fourier  components  corresponding  to  all  wavelengths  in  the  band.  When 
these  various  Fourier  components  are  added  up,  there  will  be  regions  in  the  coating  where  they 
tend  to  add  up  in  phase,  thus  producing  a  large  localized  index  modulation.  These  large 
localized  index  changes,  however,  cannot  physically  be  produced  in  the  recording  material  so  the 
intensity  recording  profile  given  in  equation  3.35  will  not,  in  reality,  always  produce  the  desired 
profile.  We  will  refer  to  this  effect  as  "index  stacking."  Although  the  amplitude  of  the 
individual  Fourier  components  making  up  the  desired  grating  structure  of  the  hologram  may  be 
small,  there  can  often  be  regions  in  the  coating  where  these  components  add  in  phase,  thus 
resulting  in  6n  values  that  are  unachievable. 


Figure  3.4  illustrates  the  concept  of  index  stacking  with  a  simple  example--a  grating 
consisting  of  two  Fourier  components  shown  in  the  top  two  curves  in  the  figure.  The  longer 
period  grating  (profile  A)  will  reflect  longer  wavelength  radiation  and  the  shorter  period  grating 
(profile  B)  reflects  shorter  wavelength  radiation.  When  the  two  index  profiles  are  added  ("A  + 
B  (Non-Chirp  Sum)")  as  shown  in  the  third  curve  in  the  figure,  the  resulting  structure  will  reflect 
both  the  long  and  short  wavelength  radiation.  Notice,  however,  the  peak-to-valley  index 
modulation  for  the  third  curve  is  equal  to  -0.35,  as  compared  to  0.20  for  the  top  two  curves 
(0. 1  separation  between  dashed  lines  in  the  figure).  When  many  different  Fourier  components 
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are  added,  index  stacking  can  present  a  real  problem. 

One  way  to  control,  and  even  eliminate,  index  stacking  is  to  judiciously  select  the  phasing 
between  the  Fourier  components  being  summed.  The  phasing  should  be  selected  to  minimize  the 
index  stacking.  In  many  cases  of  interest  we  are  only  concerned  with  the  coating  reflectivity 
profile,  R(k),  and  not  in  the  phase,  P,  of  the  reflected  radiation.  Thus  we  can  replace  r(k)  in 

equation  3.25  with  i/R (k)  exp  [  -i  P  (k)  ]  and  select  the  function  P(k)  which  minimizes  index 

stacking!  That  is,  we  must  find  the  function  P(k)  which  minimizes  the  term  in  brackets,  [  ],  in 
equation  3.25,  at  the  worst  possible  value  of  z  between  O  and  L.  Once  we  calculate  the  optimum 
P(k),  this  ideal  phase  term  can  be  used  to  find  the  desired  phase  relationship  between  the  Fourier 
components,  y  (see  equation  3.32).  This  in  turn  can  be  used  to  determine  the  location  of  the 
back  reflector  plate,  shown  in  figure  3.3,  during  the  recording  of  the  hologram.  Thus  as  the 
rotating  mirror  in  figure  3.3  is  turned  at  a  uniform  rate,  the  laser  intensity  and  back  reflector 
plate  location  are  adjusted  in  a  previously  calculated  manner.  The  changing  laser  intensity 
determines  the  amplitude  of  the  particular  Fourier  grating  components  being  recorded,  anJ  the 
back  plate  position  determines  the  relative  phase  of  that  Fourier  component. 

Although  it  is  not  obvious  from  examining  equation  3.25  that  an  ideal  P(k)  can  be  found 
which  significantly  reduces  the  effect  of  index  stacking,  some  simple  physical  reasoning  shows 
that  this  is  the  case.  This  physical  reasoning  also  can  be  used  to  directly  determine  what  the 
function  p(k)  must  be. 

The  bottommost  curve  in  figure  3.4  ("A  +  B  (Chirp  Sum)")  illustrates  the  approach  that 
can  be  used  to  minimize  index  stacking.  We  must  find  the  phasing,  P,  to  use  which  generates 
a  chirped  grating  with  constant  amplitude.  The  bottommost  curve  in  the  figure  is  a  chirped 
grating  which  contains  the  frequency  components  of  the  top  two  curves,  so  it  will  reflect  strongly 
at  the  two  desired  wavelengths.  This  curve  also  has  a  peak-to-valley  index  modulation  of  only 
0.2,  in  contrast  to  the  value  of  0.35  in  the  third  curve. 
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Our  approach  for  fabricating  "Concept  1"  holographic  filters  can  be  summarized  as 
follows: 

1 .  Calculate  the  angle  of  incidence  dependent  laser  intensity  needed  to  generate  the 
Fourier  spectrum  required  for  the  desired  R(k)  profile. 

2.  Calculate  the  Fourier  component  phasing  term,  y,  which  will  produce  a  uniform 
amplitude  chirped  grating,  as  shown  at  the  bottom  of  figure  3.4,  thus  minimizing 
the  index  stacking  problem. 

3.  Program  the  laser  output  modulator  and  the  computer  controlled  translation  stage 
holding  the  back  reflector  plate  to  respond  according  to  the  profiles  calculated  in 
steps  1  and  2. 

4.  Record  the  hologram  by  turning  the  "rotating  mirror"  in  figure  3.3  through  the 
required  range  of  angles. 

The  laser  intensity  profile  needed  for  step  1  is  given  in  equation  3.35.  The  Fourier 
phasing  term,  y,  needed  in  step  2  is  calculated  in  Appendix  3 A 

The  Fourier  component  phase  shift  term  y(0)  is  adjusted  by  precision  movement  of  the 
back  mirror  during  the  recording  process.  This  can  be  done  by  using  a  programmable  computer 
controlled  translation  stage,  such  as  Klinger  Corporation’s  MT160  Holodrive  linear  translation 
stage. 


We  introduced  the  issue  of  index  stacking  at  the  beginning  of  this  section  with  the 
example  of  a  10%  reflectivity  broad-band  (0.4  to  0.7  pm)  HC.  We  will  conclude  this  topic  by 
considering  another,  more  stressing,  coating  design.  The  desired  coating  must  have  30% 
reflectivity  from  0.367  pm  to  1.10  pm  and  0%  reflectivity  outside  of  this  band.  Since  we  do 
not  have  any  requirement  on  the  phase  of  the  reflectance,  (J(k),  this  parameter  can  be  chosen  to 
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produce  a  hologram  with  a  uniform  amplitude  index  modulation,  thus  minimizing  the  index 
stacking  problem. 

Figure  3.5  is  a  plot  of  the  index  modulation  needed  to  produce  the  wide  band  30% 
reflectivity  coating  if  we  simply  take  (5(k)  =  0.  That  is,  we  do  nothing  to  minimize  the  index 
stacking.  The  coating  has  a  total  thickness  of  20  pm  and  Figure  3.5  is  a  plot  of  the  index 
modulation  for  the  outermost  0.5  pm  of  the  structure.  The  maximum  peak-to-valley  index 
modulation  in  the  outer  portion  of  the  HC  is  nearly  equal  to  1.4  (1  +  0.4),  far  more  than  the 
0.2  maximum  which  can  be  recorded  in  DCG.  If  on  the  other  hand,  we  record  the  same  Fourier 
spectrum  as  in  figure  3.5,  but  choose  (5(k)  in  an  optimum  manner  (see  Appendix  3A),  the  peak- 
to-valley  index  modulation  can  be  reduced  to  0.2.  Figure  3.6  is  a  plot  of  the  index  modulation 
in  this  case,  for  the  outermost  0.5  pm  of  the  HC.  The  entire  profile  for  this  physically  realizable 
coating  is  a  uniform  amplitude  chirped  grating,  with  the  longest  period  portion  of  the  grating 
(which  reflects  "strongly"  (30%)  at  1.10  pm)  at  the  outermost  part  of  the  DCG  layer.  The 
period  of  the  grating  near  the  substrate  is  0. 122  pm  and  that  portion  of  the  structure  reflects  30% 
of  the  light  at  0.367  pm. 

3.3  CONCEPT  2:  DISCRETE  EXPOSURES 

3.3.1  DESCRIPTION  OF  CONCEPT  2  HOLOGRAM  FABRICATION 

Shortly  after  completing  the  analysis  for  the  "Concept  1"  hologram,  we  realized  that  a 
simpler  approach  was  available  for  fabricating  the  HCs  if  DCG  were  used.  When  a  hologram 
is  recorded  by  a  beam  at  a  single  incident  angle,  the  period  of  the  intensity  fringe  pattern,  A„, 
will  be  as  given  in  equation  3.1  (air  superstrate).  In  section  3.2  we  assumed  that  the  index  of 
refraction  modulation  of  the  developed  hologram  will  also  have  spacing  Ag.  This  holds  for  an 
ideal  recording  medium,  but  in  any  real  recording  material  nonlinearities  will  be  present  so  that 
the  index  of  refraction  profile  may  not  be  a  faithful  reproduction  of  the  material’s  exposure 
profile.  Although  DCG  is  one  of  the  best  holographic  recording  materials  discovered  to  date, 
these  non  linearities  are  present  and  prevent  a  perfect  reproduction  of  the  exposure  pattern. 
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It  turns  out  that  certain  exposure  and  development  procedures,  to  be  described  in  section 
4,  can  significantly  reduce  the  ability  of  DCG  to  faithfully  turn  a  recording  intensity  pattern  into 
a  corresponding  index  of  refraction  profile.  In  the  "Concept  2”  approach  to  matched  filter 
fabrication,  we  will  take  advantage  of  these  nonlinearities  to  fabricate  hoiugidplnc  filters  with  an 
arbitrary  reflectivity  profile.  This  technique  is  expected  to  be  simpler  than  the  "Concept  i 
approach  described  in  section  3.2. 

In  particular,  even  after  recording  a  hologram  with  a  plane  wave  at  a  single  angle  of 
incidence,  the  exposure  and  development  can  be  performed  under  conditions  that  produce  a 
chirped  grating— a  grating  with  a  period  that  changes  somewhat  through  the  thickness  of  the 
hologram.  By  carefully  controlling  the  wet  processing  of  the  exposed  DCG,  the  degree  of 
swelling  and  dehydration  of  the  developing  HC  can  be  altered.  If  the  hologram  is  fabricated  with 
a  process,  for  example,  which  causes  the  innermost  portion  of  the  coating  (near  the  substrate) 
to  swell,  the  period  of  the  index  of  refraction  modulation  will  be  greater  near  the  back  of  the 
coating  than  the  front  surface.  This  structure  will  thus  reflect  radiation  over  a  broader  bandwidth 
than  would  be  the  case  for  a  uniform  period  coating.  Examples  of  previous  applications  of  this 
property  of  DCG  will  be  given  in  section  4. 

Figure  3.7  illustrates  the  reflectivity  spectrum  of  a  uniform  period  grating  and  a  chirped 
grating.  The  spectral  response  of  the  chirped  grating  is  much  wider  than  the  uniform  period 
grating.  Also,  for  a  fixed  amplitude  modulation,  the  peak  reflectivity  of  the  chirped  structure 
will  be  less  than  the  single  period  grating. 

In  the  "Concept  2"  approach  to  holographic  filter  fabrication,  a  holographic  exposure  is 
made  at  a  number  of  discrete  incident  angles,  rather  than  over  a  continuum  of  angles.  If  the 
gratings  from  each  exposure  are  not  chirped,  the  reflectivity  of  the  HC  will  consist  of  a  series 
of  narrow  peaks— a  "comb  filter."  This  spectral  response  is  illustrated  in  figure  3.8a,  along  with 
the  continuous  profile  of  some  coating  which  we  wish  to  fabricate.  The  comb  filter  is  a  poor 
representation  of  the  desired  reflectivity  profile.  If  we  were  to  make  exposures  at  many  more 
discrete  angles,  the  resulting  hologram  would  more  closely  match  the  desired  curve.  In  the  limit 
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of  one  long  exposure  through  a  continuum  of  angles,  the  process  would  be  equivalent  to 
"Concept  I." 

Figure  3.8b  indicates  what  the  spectral  response  of  the  HC  wold  be  if  the  five  grating 
Fourier  components  which  produced  th^  comb  filter  in  figure  3.8a  were  processed  under 
conditions  which  resulted  in  chirped  gratings.  The  sum  of  the  reflectivities  from  each  of  the 
individual  chirped  gratings  can  closely  approximate  the  desired  reflectivity  profile.  If  more 
angles  are  used  for  recording  and  the  chirp  is  reduced  for  the  separate  Fourier  component,  then 
complex  reflectivity  profiles  can  be  more  accurately  fabricated.  If  too  many  separata  angles  are 
recorded,  however,  the  process  starts  to  resemble  "Concept  1"  and  the  relative  phasing  of  the 
various  gratings  will  have  to  be  controlled  in  order  to  avoid  "index  stacking."  If  the  reflectivity 
spectrum  can  be  approximated  by  only  a  few  (perhaps  less  than  10)  separately  recorded  chirped 
gratings,  then  index  stacking  will  often  not  be  a  problem. 

In  the  following  sections  we  will  analyze  in  detail  the  spectral  response  of  a  HC  fabricated 
with  the  "Concept  2"  technique.  This  analysis  will  also  help  determine  the  best  angles  and 
exposure  levels  for  recording  each  chirped  grating. 

3.3.2  CONCEPT  2  ANALYSIS 


In  order  to  determine  the  performance  of  the  "Concept  2”  HC’s  we  need  to  solve  coupled 
equations  3.16  and  3.17,  in  the  general  case  t.|B|  not  necessarily  small).  This  can  best  be 
accomplished  by  defining  a  new  variable  p: 


p(z) 


B(z) 

F(z) 


1/2 


Q(z) 


Q‘(z) 


3.38 


3-19 


where: 


K(max) 

g 

Q(z)  -  J  n,(Kg)  exp  [i  (Kg  -  2nok)z]dK„ 

..(mini 

g 

and  Q*(z)  =  complex  conjugate  of  Q(z) 

The  reflectivity  of  the  HC,  R,  is  equal  to  j  p(0">  |2. 
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If  we  take  the  derivative  with  respect  to  z  of  each  side  of  equation  3.38  and  use  equations 
3.16  and  3.17,  we  can  find  a  differential  equation  for  p(z),  which  resembles  the  Riccati 
differential  equation.  The  result  is: 


where: 


P  ' 


<T  <r 

Q  Q- 

dp 

dz 


P  '  =  -i  k  IQ  |( 


-> 

P'  + 


1)  «  4 


Q' 


-  ^ 

dz 
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The  boundary  condition  on  p(z)  is  p(L)  =  0  (L  =  coating  thickness,  z  =  0  is  at  the 
surface  of  coating  and  z  =  L  is  at  the  coating  substrate  interface).  If  the  index  of  refraction 
profile  results  from  the  superposition  of  N  different  chirped  gratings,  then  we  can  write  n:(z)  as: 


n2(z)  =  n l  +  2noY,  cos 

P"i 


K 


(  1  +  e(z))z  +  a 

(o)  P 
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In  equation  3.41  np  is  one-half  the  peak-to-valley  index  modulation  of  the  chirped 
grating  making  up  the  HC.  The  period  of  the  p^  grating  would  be  Ap(")  if  there  were  no  chirp. 
The  fractional  differential  contraction  (i.e.  the  strain)  of  the  DCG  at  coordinate  z  is  defined  as 
e(z)  (e  <  <  1),  so  that  the  local  period  of  the  p^  grating  is  equal  to  approximately  Ap""  (1  - 
e(z)).  The  term  ap  in  equation  3.41  represents  the  phase  shift  for  the  p^  chirped  grating. 

Comparing  equations  3.4  and  3.41  yields: 

..(max) 

S 

J  ft.  (K,)  exp(  i  K.z)dK, 

(min)  }  42 

S 

=  x  E  "„«p  ~  ('  ,e<z)>z  * '% 

-  r-1  A,, 

From  the  definition  of  Q,  we  find: 

,  «  \\,n  3.43 

Q(z)  =  -  V  n  exp  i  — (U  e(z))  -  2nok  z  +  i  a  ■ 

'  A  (°)  1 

P*'  [\\ 

Numerical  methods  can  now  be  used  to  solve  for  equation  3.40,  using  equation  3.43  for 
Q(z)  if  the  index  profile  can  be  expressed  as  a  sum  of  chirped  gratings  of  constant  amplitude. 
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Photo  Sensitive  (Holographic)  Medium 


Figure  3. 1  Standing  Wave  Holographic  Recording  Geometry  (Normal  Incidence) 
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—  substrate 

holographic 
material 


Figure  3.2  Geometry  of  a  Lippmann  holographic  recording  at  oblique  recording  angles. 


Laser 


Figure  3.3  Schematic  Diagram  of  Holographic  Filter  Fabrication 
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Coating  designed  for  uniform  reflectivity  of  30%  between 
the  wavelengths  of  0.367  urn  and  1.10  pm  (n0  —  1.5). 
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Holographic  Coating  with  a  Constant  Reflectance  Phase  Shift. 


Coating  designed  for  unitorm  reflectivity  of  30%  between 
the  wavelengths  of  0.367  pm  and  1.10  pm  (nQ  =  1.5). 
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Holographic  Coating  Designed  to  Minimize  "Index  Stacking." 


Figure  3.7  Reflectivity  Spectrum  of  a  Uniform  and  Chirped  Grating. 


Figure  3.8a 


Figure  3.8b 


Figure  3.8  Conventional  recording  and  development  of  a  series  of  holographic  gratings 
cannot  be  used  to  accurately  represent  a  desired  spectral  reflectivity  profile  (Figure 
3.8a).  Novel  techniques  for  recording  and  development  of  gratings  with  a  chirped 
profile  allow  the  convenient  fabrication  of  filters  with  any  desired  reflectivity 
profile  (Figure  3.8b). 
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APPENDIX  3A.  CONCEPT  1  PHASE  FACTOR  ANALYSIS 


3A.1  ANALYSIS 


In  order  to  avoid  "index  stacking,"  we  need  to  calculate  the  Fourier  component  phase 
term  y,  or  else  the  reflectance  phase  term  (J,  which  produces  a  uniform  amplitude  chirped 
grating.  The  angles  y  and  P  are  related  through  equation  3.32. 

The  square  of  the  desired  constant  amplitude  index  of  refraction  modulation  can  be  written 
as: 


n2 (z)  =  n2  +  2  no  n,  cos  g(z) 
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The  grating  vector  Kg  is  a  function  of  z  and  is  given  by  dg(z)/dz.  The  function  g(z)  is 
real  and  will  now  be  determined  in  terms  of  R(k). 

Combining  equations  3A.  1  and  3.24  yields: 
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If  we  now  multiply  each  side  of  equation  3A.3  by  exp(-i  Kg  z)  and  integrate  over  all  z 
from  to  +'»,  the  result  is: 


n,  J  exp  [ i ( g (z)  -  Kgz)]dz 
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We  will  not  show  the  details,  but  the  left  hand  side  of  equation  3A.4  can  be  approximated 
using  the  method  of  stationary  phase.  The  result  is: 
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In  equation  3A.5,  the  term  g"  is  shorthand  for  d2g/dz2  and  the  stationary  phase 
requirement  is: 

3A.6 

g'(z)  =  K 


_  dg) 

dz 


If  we  take  the  square  of  the  amplitude  of  each  side  of  equation  3A.5  and  use  equation 
3A.6,  we  find: 
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Equation  3A.7  is  our  desired  result.  It  is  a  non-linear  second  order  differential  equation 
for  the  function  g(z)  in  terms  of  known  quantities  n,,  n0,  and  R(k).  If  we  define  h(z)  as: 


2n 

h(z)  "  -uh 
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then  equation  3A.7  can  be  converted  to  a  simpler  non-linear  first  order  differential 


equation: 
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If  we  define  the  function  G(h)  in  terms  of  the  specified  reflectivity  function  as: 


G(h)  -  |  R  ( 1/h)  i 
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and  we  define  the  function  G  '(h)  as  the  inverse  function  of  G  (G  '(G(h))  =  h).  then: 
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Integrating  equation  3A.13  gives  us  our  desired  result,  the  solution  cr  g(z): 
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If  we  return  to  equation  3A.5  and  compare  the  phase  terms  on  each  side  of  the  equation. 


we  find: 
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Combining  equations  3A.  15  and  3.32,  we  can  derive  our  final  result,  an  equation  for  the 
Fourier  phase  term  y  needed  to  minimize  index  stacking: 
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y  =  g  -  z  g '  +  constant 


(The  "constant"  in  3A.  16  can  be  any  value-it  does  not  impact  the  performance  of  the  reflector) 
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In  order  to  properly  control  the  position  of  the  back  reflector  plate  during  the  recording 
of  the  HC,  we  must  actually  know  y  as  a  function  of  0,  the  recording  beam’s  angle  of  incidence 
on  the  hologram  during  the  angle  sweep.  This  is  easily  found  by  the  following  process: 

1.  Use  equation  3A.14  to  solve  for  g(z) 

2.  Use  the  relationship 

g'(z)  =  Kg(0)  =  (nj  -  n^sin20)1/2 

r 

to  solve  for  z  =  z  (0) 

3.  Use  z(0)  in  the  right  hand  side  of  equation  3A.16  to  und  y(0). 

3A.2  SAMPLE  CALCULATION 


We  will  now  apply  the  results  of  the  last  section  to  the  simple  case  of  a  HC  with  a 
constant  reflectivity  R(,  between  the  wavelengths  A.min  and  kmax.  Since  R0  is  a  constant,  equation 
3A.  10  gives  the  result: 
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The  constant  c,  is  easily  found  by  evaluating  equation  3A.  18  at  the  point  z  =  0.  We  thus 

find: 
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If  we  replace  the  left  hand  side  of  equation  3A.19  with  K„(0)  and  solve  for  z,  the  result 
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Since  y  =  g  -  zg'  (equation  3A.16),  we  can  use  equations  3A.20,  3A.21,  and  g'  =  Kg 
to  find  y. 
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The  function  K  (0)  is  given  in  equation  3.28. 
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4.  PROPERTIES  OF  HOLOGRAPHIC  RECORDING  MATERIALS 


4.1  BACKGROUND 


Dichromated  gelatin  (DCG)  is  probably  the  best  material  for  fabricating  high  quality 
holograms.  It  offers  one  of  the  largest  available  refractive  index  modulations  (  —  0.1),  has 
excellent  resolution  (>  10,000  lines/mm),  has  low  scattering,  and  requires  moderate  exposures 
(50  -  300  mJ/cm2  in  the  wavelength  range  350  -  520  nm).  It  can  be  deposited  in  layers  up  to 
50  microns  thick,  which  is  equivalent  to  over  300  periods  for  recording  with  radiation  at  0.5  pm. 
The  ability  of  DCG  to  provide  an  index  modulation  of  0. 1  in  layers  up  to  50  pm  thick  is  an 
important  feature  for  our  intenced  application.  Figure  4.1  is  a  table  comparing  the  properties 
of  four  different  materials  which  have  been  used  for  holographic  recording.  The  high  resolution 
of  dichromated  gelatin  combined  with  its  large  index  modulation  make  it  the  primary  material 
for  the  Phase  II  work. 

Another  desirable  property  of  DCG  is  its  high  optical  damage  threshold.  JBS 
Technologies,  Inc.  received  a  $103,000  subcontract  from  SAIC,  as  part  of  an  SDIO  coating 
development  program,  to  investigate  the  potential  of  HCs  for  near  IR  lasers. r>  Key  findings  of 
this  work  were  that  damage  thresholds  of  over  1  MW/cm2  cw  were  achievable  for  HCs  deposited 
on  Moly  and  that  reflectivities  of  99.5%  were  possible  at  a  wavelength  of  1.3  pm.'11 

4.2  THE  DCG  RECORDING  PROCESS 


The  unsensitized  recording  material  is  gelatin,  which  consists  of  the  fibrous  protein 
collagen;  a  material,  normally  obtained  from  animal  skins  and  bones.  Gelatin  has  the  property 
that  it  can  be  cross-linked  by  certain  polyvalent  metallic  ions  when  exposed  to  light.  The 
chromium  ion  (CiJ+)  is  the  most  commonly  used  metal  ion  for  this  purpose.  It  is  obtained  from 


nSAIC  contract  DAAH01-89-R-0065,  D.O.  #1  to  the  Army  (supported  by  SDIO). 
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a  reduced  inorganic  dichromate  which  reacts  at  the  free  carboxyl  group  sites  in  gelatin,  provided 
by  the  constituent  di-basic  aspartic  and  glutamic  amino  acids.  The  photo-reduction  of  the 
dichromate  ion  and  its  effect  on  the  crosslinking  of  gelatin  through  the  carboxyl  functional  groups 
is  illustrated  in  figure  4.2® 

This  crosslinking  significantly  increases  the  molecular  weight  and  reduces  the  solubility 
of  the  material  in  warm  water.  Thus,  portions  of  the  DCG  exposed  to  radiation  in  the  range 
3500  to  5200A  will  become  less  water  soluble  (hardened)  and  this  hardening  will  be  proportional 
to  the  exposure  level.  When  the  material  is  then  washed  and  subsequently  dried  in  a  dehydrating 
agent,  such  as  isopropanol,  the  gelatin  layer  will  have  a  varying  refractive  index  structure 
recorded  in  it.  This  index  of  refraction  pattern  corresponds  to  the  exposure  levels  which  various 
parts  of  the  DCG  were  subjected  to  during  the  holographic  recording  process.  Figure  4.3 
summarizes  the  DCG  recording  process.® 

Of  particular  interest  in  the  current  work  is  how  the  index  of  refraction  modulation  pattern 
can  be  modified  from  the  sinusoidal  standing  wave  intensity  pattern  of  the  recording  radiation. 
In  other  words,  how  a  chirped  or  tapered  grating  can  be  created  from  the  purely  sinusoidal 
exposure  pattern.  There  are  four  main  stages  in  the  DCG  recording  procedures:  (1) 
preprocessing;  (2)  sensitization;  (3)  exposure  and  (4)  development  or  processing  (water  washed 
followed  by  dehydration).  It  is  primarily  the  procedure  followed  in  the  final  processing  stage 
which  can  be  used  to  control  grating  chirp.  We  will  discuss  this  step  in  detail  in  the  following 
section.  A  description  of  the  preprocessing,  sensitization,  and  exposure  steps  will  not  be  given 
here,  but  can  be  found  in  references  by  Cullen,®  Kubota, (3)  Sjolinder,®  and  Georgekutty.® 

4,3  DCG  PROCESSING  PROCEDURE  TO  CONTROL  GRATING  CHIRP121 

The  development  or  processing  step  begins  with  a  water  wash  to  remove  the  unreacted 
ammonium  dichromate  after  the  exposure  of  the  DCG  Increasing  the  water  temperature  and  pH 
will  increase  the  swelling  of  the  gelatin  layer  which  increases  the  index  of  refraction  modulation. 
Gelatin  molecules  are  transferred  and  rearranged  in  the  water  swollen  layer  to  regions  of  higher 
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exposure,  as  illustrated  in  figure  4.3.  Since  the  absorbed  water  produces  swelling  in  the  DCG, 
the  period  of  the  index  modulation  is  greater  than  that  of  the  sinusoidal  intensity  pattern  created 
during  the  exposure  step.  The  hologram,  at  this  stage  of  the  development  process,  would  thus 
most  strongly  reflect  radiation  at  a  wavelength  slightly  longer  than  the  recording  wavelength. 
The  procedure  followed  in  the  next  step  of  the  processing,  dehydration,  can  be  used  to  freeze 
in  this  swelling,  reverse  it,  or  to  retain  the  swelling  in  certain  parts  of  the  DCG  layer  and  reverse 
it  in  other  parts.  It  is  this  third  possibility  that  produces  chirped  gratings  from  a  purely 
sinusoidal  recording. 

Isopropanol  is  considered  to  be  the  best  dehydrating  agent  for  DCG.  When  the  sample 
is  immersed  in  the  isopropanol,  the  water  in  the  swollen  gelatin  is  rapidly  removed,  thus  freezing 
the  DCG  molecular  structure.  If  100%  isopropanol  is  used  for  the  dehydration  bath  immediately 
after  the  water  bath,  the  dehydration  will  occur  so  rapidly  at  the  surface  that  some  water  will  be 
trapped  deeper  in  the  DCG  layer-closer  to  the  substrate.  This  water  is  unable  to  escape  through 
the  dry  hardened  surface  and  the  DCG  layers  closer  to  the  substrate  will  remain  permanently 
swollen.  Thus  the  coating  will  naturally  be  chirped,  with  the  longer  period  closer  to  the 
substrate. 

An  unchirped,  faithful  reproduction  of  the  sinusoidal  recording  intensity  pattern  can  be 
generated,  on  the  other  hand,  if  the  rate  of  dehydration  is  carefully  controlled.  A  slower  rate 
of  dehydration  results  if  graduated  water/alcohol  mixtures  are  used.  The  sample  is  taken  from 
the  wa-er  bath  and  initially  immersed  in  a  90/10  solution  of  water/alcohol.  It  can  then  be  added 
to  a  50-  50  mixture  of  water  and  alcohol,  followed  by  a  10/90  solution,  and  finally  a  bath  of  pure 
isoprcyanol.  This  allows  sufficient  time  for  a  more  uniform  dehydration  of  the  developing  DCG 
and  a  r'nusoidal,  rather  than  chirped,  index  of  refraction  modulation  will  be  achieved.  The  level 
of  chir.iing  can  thus  be  carefully  controlled  by  controlling  the  rate  of  dehydration.  Table  4. 1 
summa  izes  typical  processing  procedures  for  making  narrow  bandwidth  (no  chirp)  and  wide 
bandwidth  (chirped)  coatings. 
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Table  4.1  Typical  Processing  Steps  for  Sinusoidal  and  Chirped  Coatings 


Narrow  Bandwidth  Process  Wide  Bandwidth  Process 

•  Prepare  water/IPA  mixtures,  90/10,  50/50,  Soak  in  100%  IPA  (20-25°C) 

10/90,  0/100%  and  soak  plate  in  each  with  agitation  until  maximum 

(starting  with  the  90/10%)  for  4  min  efficiency  is  obtained, 

(20-25 °C)  with  agitation.  typically  5-8  min. 

•  Processed  plates  are  transferred  to  a  controlled  RH  environment  of  approx.  40%  to 
stabilize.  10-20  hours. 

•  To  protect  the  hologram  against  humidity  changes  a  cover  glass  can  be  adhered  using 
suitable  transparent  cyanoacrylate  or  epoxy  resins. 


It  has  been  demonstrated  that  the  shape  of  the  reflectivity  profile,  as  well  as  its  width,  can 
be  controlled  through  careful  processing  and  variations  in  the  gelatin  formulation.  Figure  4.4 
illustrates  how  Cullen(2)  was  able  to  fabricate  a  coating  with  a  pointed  reflectivity  peak  or  a  flat 
top,  depending  on  the  processing  used. 

The  author  of  this  report  has  also  needed  to  fabricate  chirped  HCs  for  the  previously 
referenced  SDIO  program.0*  The  holographic  coatings  for  that  program  were  needed  for  a 
chemical  laser  (HF  overtone)  with  a  central  wavelength  of  1.35  pm  and  a  bandwidth  of  a  little 
over  0.1  pm.  By  following  the  development  procedures  discussed  earlier,  a  broad-band  high 
reflectivity  coating  with  a  flat  top  profile  was  produced.  Figure  4.5  is  a  plot  of  the  transmission 
and  optical  density  of  one  of  the  fabricated  coatings. 

4.4  SEQUENTIAL  RECORDINGS  IN  DCG 


Sequential  recording  of  different  gratings  in  a  single  DCG  hologram  have  been  reported 
by  numerous  groups  such  as  Magarifios  and  Coleman/6’  Rallison/7*  and  Physical  Optics 
Corporation/8’  In  the  work  reported  by  Rallison,  five  overlapping  gratings,  each  with  near  100% 
diffraction  efficiency,  were  recorded  in  DCG  20  pm  thick.  In  another  device  fabricated  by 
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Rallison,  nine  sequential  exposures  were  made  on  a  single  hologram.  The  nine  exposures  were 
of  lenses  focusing  a  plane  input  beam  to  different  locations.  When  the  developed  hologram  was 
subsequently  illuminated  by  a  plane  wave,  the  wave  emerging  from  the  back  of  the  hologram 
focused  down  to  nine  diffraction  limited  spots.  The  nine  separate  gratings  were  recorded  so  that 
the  energy  in  the  nine  output  spots,  making  up  the  3  x  3  array,  was  evenly  distributed.  Rallison 
has  stated  that  the  maximum  number  of  separate  high  efficiency  recordings  that  can  be  made  in 
a  single  hologram  is  about  equal  to  twice  the  index  modulation  times  the  thickness  in  microns. 
Thus  a  50  micron  thick  layer  of  DCG  should  be  sufficient  for  recording  ten  separate  high 
efficiency  chirped  gratings,  which  when  pieced  together  can  form  the  desired  reflectivity 
spectrum.  Rallison’s  estimate  is  somewhat  qualitative, (,)  and  a  much  more  precise  prediction 
can  be  made  with  the  computer  codes  written  for  this  contract  and  discussed  in  the  next  section. 
These  codes  are  based  on  the  analysis  presented  in  section  3  of  this  report. 

Although  many  groups  have  fabricated  chirped  gratings  in  DCG  and  a  number  of  groups 
have  fabricated  HCs  with  multiple  gratings,  none  have  combined  these  two  features  in  order  to 
make  a  coating  with  some  prescribed  reflectivity  spectrum.  The  separately  recorded  coatings 
have  nearly  always  been  for  wavelengths  well  separated,  so  that  interaction  effects  could  be 
ignored.  In  our  case,  we  need  to  combine  the  spectral  profiles  of  each  recorded  grating  to  form 
a  variety  of  possible  reflectivity  profiles. 


nFor  example,  what  exactly  is  meant  by  "high  efficiency"  -  70%,  90%,  99%  diffraction 
efficiency? 
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Figure  4.1.  Comparison  of  Four  Materials  Used  in  Holographic  Recording. 
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Figure  4.2  Photo-reduction  of  the  dichromate  ion  and  its  effect  on  the  crosslinking  of  gelatin 
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Figure  4.3  Generalized  DCG  recording  process'2' 
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Figure  5.  Measured  Transmission  and  Optical  Density  of  a  Lippman  Mirror  on  Glass 
Substrate. 
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5.  COMPUTER  MODELING  OF  OPTICAL  COATING  PERFORMANCES 


A  computer  program  has  been  written  to  find  numerical  solutions  to  the  coupled  mode 
equations.  The  program  uses  a  subroutine  based  on  the  fifth  order  Runge-Kutta  method  w'ith 
adaptive  step  size.  Appendix  5A  is  a  listing  for  this  program. 

The  program  was  used  to  design  three  different  filters: 

1 .  A  filter  with  a  reflectivity  profile  matching  the  UV  spectral  signature  (3000A  to 
4000A)  of  a  typical  solid  propellant  booster,  as  shown  ir.  figure  2.9. 

2.  A  filter  with  a  reflectivity  profile  given  by  the  optimum  matched  filter 
characteristics  used  to  detect  a  solid  fueled  booster  against  an  Earth  background 
(mid-latitude  summer  profile  illustrated  in  figure  2.8). 

3.  High  uniform  reflectivity  (ideally  equal  to  1.0)  over  the  waveband  0.4  to  0.8  qm 
and  zero  elsewhere. 

The  first  two  filters  listed  above  would  be  used  in  a  reflectivity  mode  -  they  would  be 
placed  in  the  camera  beam  train  so  that  light  reflected  from  them  would  be  focused  on  the 
detector  used  to  detect  the  booster  target.  They  would  preferen*  illy  reflect  radiation  at  the 
stronger  spectral  lines  of  the  target  and  reflect  a  smaller  fraction  of  Earth  background  radiation 
onto  the  detectors. 

5.1  UV  PLUME  SPECTRAL  PROFILE 


Figure  5.1a  illustrates  the  UV  spectral  signature  of  a  solid  booster,  while  figures  5.1b, 
c,  d,  and  e  show  the  calculated  reflectivity  profiles  of  four  different  HCs  fabricated  in  DCG  with 
slightly  different  parameters.  All  four  reflectivity  profiles  are  reasonably  accurate  representations 
of  the  plume  signature,  but  the  profiles  in  figure  5.  Id  and  5.  le  are  the  best  representations.  The 
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design  parameters  for  the  four  different  reflectors  are  given  in  figure  5.2.  All  require  very  small 
chirp  parameters  ( 1 .5%),  so  the  slow,  well  controlled  dehydration  procedure  described  in  section 
4.3  would  be  used  to  fabricate  the  holographic  reflectors  in  DCG.  The  exposure  time  for  each 
of  the  four  separate  recordings  will  be  proportional  to  the  desired  reflectivity.  We  have  found 
that  slightly  better  representations  of  the  four  peaks  shown  in  figure  5.1a  can  be  realized  if  the 
side  lobes,  that  are  normally  present  in  fairly  narrow  band  reflectors  (1.5%  chirp),  can  be 
suppressed.  This  can  be  accomplished  by  recording  the  four  gratings  in  a  manner  that  produces 
a  tapered  index  modulation.  The  index  modulation  nt  (i  =  1  to  4)  used  to  arrive  at  the  four 
profiles  shown  in  figure  5.1  was  quadraecally  tapered  and  of  the  form: 


(max) 

n 
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( max ) 

n. 


(min)\ 
"■  ) 


z 
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where: 


I 

n  (max)  _ 


n  (mm)  _ 


DCG  thickness 

maximum  index  modulation  for  the  i—  recording  (occurs  at  z  =  0 
and  L,  the  hologram  surface  and  DCG/substrate  interface) 
minimum  index  modulation  for  the  i—  recording  (occurs  at  z  = 
J2) 


The  index  modulation  taper  can  be  adjt  .ted  through  the  selection  of  recording  wavelength 
and  preparation  of  the  DCG  (references  2  and  4  from  section  4  of  this  report).  As  mentioned 
above,  however,  the  curses  plotted  in  figure  5. 1  are  not  too  sensitive  to  this  parameter,  so  it  will 
not  have  to  be  controlled  too  carefully  during  HC  fabrication.  The  values  of  n,(ITUn)  and  n,(max) 
are  given  in  figure  5.2. 

The  HC  thickness  needed  to  achieve  the  four  profiles  shown  in  figure  5.1  is  only  10  pm 
and  the  parameter  a,  -  2n()A1<°>  for  the  Hnir  recordings  is  also  given  in  figure  5.2 


5.2  VISIBLE  PLUME  SPECTRAL  PROFILE 

Figure  5.3a  illustrates  the  optimum  matched  filter  profile  for  the  detection  of  a  solid 
booster  against  an  Earth  background,  as  calculated  in  section  2.1.5.  Figures  5.3b,  5.4a,  and 
5.4b  are  plots  of  the  reflectivity  spectrum  of  several  HCs  designed  to  match  the  optimum  filter. 
The  coating  illustrated  in  figure  5.3b  is  the  best  representation  of  the  desired  spectral  profile. 
The  parameters  used  to  calculate  these  profiles  are  given  in  figure  5.5.  All  three  HCs  require 
a  coating  thickness  of  22  pm,  a  small  chirp  (1  %),  and  only  two  separate  recordings  (two  chirped 
coatings  recorded  in  the  DCG)  to  produce  the  indicated  spectral  profiles. 

5.3  HIGH  REFLECTIVITV  COATING  FOR  THE  VISIBLE 

Finally,  a  much  more  difficult  spectral  profile  to  produce  holographically  was  the  case 
of  a  reflectivity  of  1.0  from  0.4  pm  to  0.3  pm,  and  zero  elsewhere.  We  attempted  to  design  a 
HC  with  this  profile  using  two  and  then  four  separate  recordings.  The  total  index  modulation 
of  the  DCG,  consisting  of  the  ium  of  the  modulation  from  each  separate  recording,  was 
constrained  to  have  a  peak-to-vallev  value  no  greater  than  0.2  (see  section  4.1).  The  maximum 
allowable  fractional  chirp  was  constrained  to  be  no  greater  than  15%,  a  conservative  assumption 
for  DCG. 

Figure  5.6  shows  the  reflectivity  profiles  of  HCs  recorded  with  two  different  chirped 
periods.  The  input  parameters  are  summarized  in  figure  5.7.  Figure  5.6a  used  a  5%  chirp, 
which  is  seen  to  be  totally  inadequate  for  spanning  the  wide  spectral  band  of  interest. 
Significantly  better  results  are  achieved  by  combining  a  10%  chirp  with  a  smaller  k2  value  (0.75 
pm  -  0.60  pm),  as  illustrated  in  figures  5.6b  and  5.6c.  Surprisingly,  increasing  the  chirp  from 
10%  to  15%,  while  keeping  other  parameters  fixed  (figure  5.6c  -  5.6d)  did  not  improve  the 
profile.  The  performance  dropped  somewhat  in  figures  5.6e  and  5.6f,  as  compared  to  figure 
5.6d,  because  the  coating  thickness  was  reduced  from  50  pm  to  40  pm.  In  figures  5.6g,  h,  and 
i  the  DCG  thickness  is  increased  to  50  pm  and  the  average  period  of  the  shorter  period  grating 
is  increased  somewhat.  This  fills  in  the  reflectivity  drop  centered  at  0.57  pm  in  figures  5.6b 


5-3 


through  5.6f,  promoting  a  higher  more  uniform  reflectivity  profile.  It  does,  however,  result  in 
a  short  wavelength  spectral  cutoff  at  around  0.45  to  0.5  pm,  rather  than  the  desired  value  of  0.4 
pm. 


Figure  5.8a  through  5.8j  illustrate  the  reflectivity  profile  when  four  separate  recordings 
are  made  in  the  DCG  to  generate  four  separate  chirped  grating  profiles.  The  coating  input 
parameters  for  these  ten  cases  are  summarized  in  figure  5.9.  These  reflectivity  profiles  are  not 
as  close  to  our  design  goal  as  was  the  case  for  the  two  exposure  profiles  given  in  figure  5.6. 
This  is  primarily  due  to  the  fact  that  a  smaller  index  modulation  per  recording  (n(  factor)  had  to 
be  used  for  the  four  recording  case,  in  order  to  avoid  "index-stacking"  (ensure  that  the  peak-to- 
valley  modulation  was  never  more  than  0.21. 

The  15%  maximum  chirp  profile  and  50  pm  DCG  thickness  assumed  in  generating  the 
best  eflectivity  profiles  for  the  N  =  2  (two  exposure)  case  may  be  a  bit  conservative.  It  may 
be  possible  to  produce  a  slightly  larger  chirp  (20%)  on  thicker  DCG  layers  ( -  70  pm),  in  which 
case  a  high  reflectivity  coating  from  0.4  to  0.8  pm  could  be  produced,  based  on  the  trends  shown 
in  figure  5.6. 

5-4  PUBLICATION  OF  CONTRACT  CONCLUSIONS 

A  summary  of  the  results  of  this  report,  including  the  "plume  filter"  examples  in  section 
5. 1  and  5.2,  have  been  submitted  for  publication  and  presentation  at  the  1991  SPIE  International 
Symposium  on  Optical  Applied  Science  and  Engineering  (21-26  July  1991,  San  Diego,  CA; 
Conference  No.  15).  A  draft  copy  of  the  paper  is  given  in  Appendix  5B. 
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OPTIMUM  MATCHED  FILTER  REFLECTIVITY  AS  A  FUNCTION  OF 
WAVELENGTH  FOR  THE  DETECTION  OF  A  SOLID  FUELED  BOOSTER 
AGAINST  AN  EARTH  BACKGROUND  USING  A  SILICON  DETECTOR 


5.3a 


5.3b 

I  igure  5.3  Optimum  Matched  Filter  Reflectivity,  as  Calculated  for  a  Solid  Fueled 

Against  an  Earth  Background  (5.3a)  and  the  Reflectivity  Profile  of  an  HC  (5.3  -n 
with  Parameters  Given  in  Figure  5.5. 


[  igure  5.4 


Holographic  Matched  Filters  for  the  Detection  ot  a  Solid  Fueled  Booster  Again  v. 
an  Earth  Background,  with  Coating  Parameters  Given  in  Figure  5.5.  These 
Filters  do  not  Reproduce  the  Optimum  Profile  (Figure  5.3a)  Quite  as  well  as  the 
HC  Shown  in  Figure  5.3b. 
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Figure  5.9  HC  parameters  used  to  generate  the  reflectivity  plots  in  figure  5. 8  (four  chirped  coatings;  maximum  allowable  peak-to-valley  index  of  refraction 
modulation  =  0.2). 


COUPLED  MODE  EQUATIONS 


program  CaseTwo 
c 

c  This  is  a  program  for  calculating  the  reflectivity 

c  profile  of  multiple  chirped  gratings, 

c  The  program  uses 

c  Runge-Kutta  driver  with  adaptive  stepsize  control, 

c  Integrate  the  NVAR  starting  values  YSTART  from  XI  to  X2 

c  with  accuracy  EPS,  storing  intermediate  results  in  the 

c  OUTPUT  file.  H 1  should  be  set  as  a  guessed  first 

c  stepsize.  HMIN  as  the  minimum  allowed  stepsize  (can  be  0). 

c  YSTART  is  replaced  by  values  at  the  end  of  the 

c  integration  interval.  DERIVS  is 

c  the  subroutine  for  calculating  the  right-hand 

c  side  derivative,  while  RKQC  is  the  name  of  the  stepper 

c  routine  to  be  used, 

c 

complex  ystart(  1 7),y(  17),dydx(  17) 
real  yscal(  17) 
real  ynorm(2) 
real  xl,wvl,wvk,thk 
real  zeron 
real  pi 
complex  ci 
c 

common/mainp/x  1  ,wvl,wvk,thk 
common/allcom/zeron,pi,ci 
c 

c  . parameters . 

two=2.0 

zero=0.0 

tiny=l.e-30 

call  params(nwvl,wvlmin,wvlmax,delwvl) 

c  . parameters . 

c 

open(2,file='case2.d  1  ’) 
c 

hmin=l.e-7 

eps=l.e-3 


c  write(*,*)'input  cutoff  factor  (0.5  ?)’ 

c  read(*,*)ctff 

ctff=0.5 
c 

c  . start  of  the  big  loop . 


do  90  iwvl=l,nwvl 
wvl=wvlmin+delwvl*float(iwvl- 1 ) 
wvk=2.*pi*zeron/wvl 
thk=wvk*ctff 

c  maxstp=max.  number  of  points  from  xl=L  to  x2=0 

maxstp=  10000 
c  L=10  micron 

xl=10. 
x2=0. 
nvar=2 
dxsav=0.01 
hi  =0.01 
x=xl 
h=hl 

ystart(l)=(l.,0.) 
ystart(2)=(0.,0.) 
do  11  i=l,nvar 

y(i)=ystart(i) 

1 1  continue 
xsav=x-dxsav*two 
do  16  nstp=l,maxstp 

call  derivs(x,y,dydx) 
c  write(*,*)nstp 

do  12  i=l,nvar 

yscal(i)=cabs(y(i))+cabs(h*dydx(i))+tiny 

12  continue 

if((x+h-x2)*(x+h-x  1  ).gt.zero)  h=x2-x 

call  rkqc(y,dydx,nvar,x,h,eps,yscal,hdid,hnext,derivs) 

if(abs(x-xsav).ge.dxsav)  then 

xsav=x 

endif 

if((x-x2)*(x2-xl).ge.zero)  then 
do  14  i=l,nvar 

ystart(i)=y(i) 


14 

continue 
do  100  i=l,nvar 

ynorm(i)=cabs(y(i)) 

100 

continue 

yt  1  =ynorm(  1  )*ynorm(  1 ) 
yt2=ynorm(2)*ynorm(2) 
R=yt2/ytl 

write(2,5000)wvl,R 

format(lx,g20.10,' 


5000 


,g20.10) 


goto  90 

endif 

if(abs(hnext).lt.hmin) 

*  pause  'Stepsize  smaller  than  minimum.' 
h=hnext 

16  continue 

pause  Too  many  steps.' 

90  continue 

stop 
end 
c 
c 
c 

subroutine  rkqc(y,dydx,n,x,htry,eps,yscal,hdid,hnext,derivs) 
c  Fifth-order  Runge-Kutta  step  with  monitoring  of  local 

c  truncation  error  to  ensure  accuracy  and  adjust  stepsize. 

c  Input  are  the  dependent  variable  vector  Y  of  length  N 

c  and  its  derivative  DYDX  at  the  starting  value  of  the 

c  independent  variable  X.  Also  input  are  the  stepsize  to  be 

c  attempted  HTRY,  the  required  accuracy  HNEXT  is  the 

c  estimated  next  stepsize.  DERIVS  is  the  user-supplied 

c  subroutine  that  computes  the  right-hand  side  derivatives, 

parameter  (pgrow=-0.2,pshmk=-0.25,fcor=l./15., 

*  one=l.,safety=0.9,errcon=6.e-4) 

c  The  value  ERRCON  equals  (4/SAFETY)* *(l/PGROW) 

external  derivs 

complex  y(  1 7),dydx(  1 7),ytemp(  1 7),ysav(  1 7),dysav(  1 7) 
real  yscal(17) 
xsav=x 
do  1 1  i=l,n 

ysav(i)=y(i) 

dysav(i)=dydx(i) 

1 1  continue 

h=htry 

1  hh=0.5*h 

call  rk4(ysav,dysav,n,xsav,hh,ytemp, derivs) 
x=xsav+hh 

call  derivs(x,ytemp,dydx) 

call  rk4(ytemp,dydx,n,x,hh,y, derivs) 
x=xsav+h 

if(x.eq.xsav)pause  'stepsize  not  significant  in  RKQC' 
call  rk4(ysav,dysav,n,xsav,h,ytemp, derivs) 
errmax=0. 


1 

1 

1 

do  12  i=l,n 

ytemp(i)=y(i)-ytemp(i) 

errmax=max(errmax,cabs(ytemp(i)/yscal(i))) 

1 

1 

12 

continue 

errmax=errmax/eps 

if(errmax.gt.one)then 

h=safety*h*(errmax**pshmk) 
goto  1 
else 

hdid=h 

1 

1 

1 

if(errmax.gt.errcon)then 

hnext=safety*h*'(errmax**pgrow) 

else 

hnext=4.*h 

endif 

endif 

do  13  i=l,n 

y(i)=y(i)+ytemp(i)*fcor 

1 

13 

continue 

■ 

19 

return 

1 

m 

c 

c 

c 

end 

i 

subroutine  rk4(y,dydx,n,x,h,yout,derivs) 

c 

Given  values  for  N  variables  Y  and  their  derivatives 

■ 

c 

DYDX  known  at  X,  use  the  fourth-order  Runge-Kutta  method 

1 

c 

to  advance  the  solution  over  an  interval  H  and  return  the 

c 

incremented  variables  as  YOUT,  which  need  not  be  a  distinct 

■ 

c 

array  from  Y.  The  user  supplies  the  subroutine 

■ 

c 

DERIVS(X,Y,DYDX)  which  returns  derivatives  DYDX  at  X. 
external  derivs 

i 

i 

■ 

complex  y(  17),dydx(  1 7),yout(  1 7),yt(  1 7),dyt(  1 7^  d vm(  17) 

hh=h*0.5 

h6=h/6. 

xh=x+hh 

do  11  i=l,n 

yt(i)=y(i)+hh*dydx(i) 

1 

i 

11 

continue 

call  derivs(xh,yt,dyt) 
do  12  i=l,n 

yt(i)=y(i)+hh*dyt(i) 

i 

i 

12 

continue 

call  derivs(xh,yt,dym) 
do  13  i=l,n 

yt(i)=y(i)+h*dym(i) 

dym(i)=dyt(i)+dym(i) 

continue 

call  derivs(x+h,yt,dyt) 
do  14  i=l,n 

yout(i)=y(i)+h6*(dydx(i)+dyt(i)+2.*dym(i)) 

continue 

return 

end 


subroutine  params(nwvl,wvlmin,wvlmax,delwvl) 
real  alpha 

real  phi(  1 5),grtk(  15) 
integer  ngrt 
real  zeron 
real  pi 
complex  ci 

common/prmtrs/alpha, phi, grtk, ngrt 
common/allcom/zeron,pi,ci 

. parameters  . 

ci=(0.,l.) 
pi=3. 14159 
zeron=n_0  : 
zeron=1.5 

alpha=absorption  constant  (  1 /micron)  : 
alpha=0. 

ngrt=number  of  gratings  (N)  : 
ngrt=4 

grating  phase  shift  : 

phi(l)=0. 

phi(2)=pi 

phi(3)=0. 

phi(4)=0. 

grating  wave  length  (micron)  grtk=4  n_0  pi/wvl: 
grtk(l)=4.*zeron*pi/0.315 
grtk(2)=4.*zeron*pi/0.3365 
grtk(3)=4.*zeron*pi/0.3580 


c 

c 

c 


grtk(4)=4.*zeron*pi/0.39 
. parameters . 


wvl=lambda  (wavelength)  micron 
nwvl=60 
wvlmin=0.3 
wvlmax=0.4 

delwvl=(wvlmax-wvlmin)/float(nwvl- 1 ) 
return 
end 
c 
c 

subroutine  derivs(x,y,dydx) 
complex  y(17),dydx(17) 
complex  g,gcc 
real  deln(15),epsilon(15) 
real  alpha 

real  phi(15),grtk(15) 
integer  ngrt 
real  zeron 
real  pi 
complex  ci 
real  xl,wvl,wvk,thk 
c 

common/mainp/x  1  ,wvl,wvk,thk 
common/prmtrs/alpha, phi, grtk, ngrt 
common/allcom/zeron,pi,ci 


c  . parameters . 

c  deln=n_l  (tapered  if  a  function  of  x) : 

deln(  1  )=0.009+0.008*x*x/(x  1  *x  1  )-0.008*x/x  1 
deln(2)=0.0065+0.004*x*x/(x  1  *x  1  )-0.004*x/x  1 
deln(3)=0.0085+0.004*x*x/(x  1  *x  1  )-0.004*x/x  1 
deln(4)=0.01 8+0.007  *x*x/(x  1  *x  1  )-0.007*x/x  1 
c  epsilon=chirp  factor 


epsilon(  1  )=  1  .-0.0 1 5*x/x  1 
epsilon(2)=  1  .-0.0*x/x  1 
epsilon(3)=  1  .-0.0*x/x  1 
epsilon(4)=  1  .+0.0 1 5*x/x  1 

g=(0.,0.) 

do  10  i=l,ngrt 
dk=2.*wvk-grtk(i) 
if(abs(dk).ge.thk)  goto  10 


delk=2.*wvk-epsilon(i)*grtk(i) 
g=coefficient  in  the  coupled  mode  equation 
g=g+(ci*pi*deln(i)/wvl)*cexp(ci*(delk*x-phi(i))) 
continue 
gcc=conjg(g) 

. parameters . 

dydx(  1  )=g*y(2)-0.5  *alpha*y(  1 ) 
dydx(2)=gcc*y(  1  )-0.5*alpha*y(2) 
return 
end 
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ABSTRACT 

We  analyze  the  reflection  spectrum  of  multiple  chirped  gratings  using  coupled- 
inode  theory.  We  also  design  filters  with  desired  reflection  spectrum  by  adjusting 
the  parameters  of  multiple  chirped  gratings. 


I 


1.  Introduction 


It  is  important  to  design  filters  with  desired  reflection  spectrum.  These  filters 
require  certain  variation  of  index  of  refraction  in  a  medium.  Holography  is  one  of 
the  most  practical  methods  to  generate  any  desired  index  variation  in  a  medium. 
In  this  paper,  we  report  the  design  of  filters  with  multiple  chirped  holographic 
gratings.  We  first  analyze  the  reflection  spectrum  of  multiple  chirped  gratings. 
Then,  we  adjust  different  parameters  to  design  filters  with  desired  reflection  spec¬ 
trum. 


2.  Coupled-mode  Analysis 


Referring  to  Fig.  1,  a  beam  of  light  is  incident  upon  a  medium  with  thickness 
L.  The  transmittance  and  reflectivity  depend  on  both  the  index  of  refraction  of 
the  medium  and  the  wavelength  of  light.  The  reflection  spectrum  can  be  obtained 
by  calculating  the  reflectivity  of  plane  waves  with  different  wavelength. 


Consider  a  medium  where  the  index  variation  is  generated  by  a  set  of  N 
holographic  gratings.  The  total  index  of  refraction  can  be  written  as 

•v 


n(c)  =  n0  +  Y"  rip(c) cos[^-(l  +  e(c))c  +  op  +  o], 
P=I  p 


r 


where 


np(s)  =  4°>  +  np]J  +  np2)~[J 


(2) 


is  the  amplitude  of  each  grating  which  is  tapered  during  manufacturing,  f  is  the 
chirp  factor,  L  is  the  thickness  of  the  medium,  A;)  is  the  grating  spacing  of  the 


•2 


p-th  grating,  4>p  and  a  are  the  phase  and  absorption  constant  respectively.  A 
plane  wave  incident  upon  the  medium  will  be  reflected.  The  total  electric  field 
inside  the  medium  is 


E  =  A{z)e,{ijt-k:)  +  B(z)e'{ujt+k:). 


(3) 


where  A(z)  and  B(z)  represent  the  amplitudes  of  the  forward  and  backward 
propagating  plane  waves  respectively.  Under  the  approximation  of  slow  varying 
envelopes,  the  coupled-mode  equations  [1]  can  be  written  as 


d_A  =  iy^Memzk-K,u-.,\D  _  i„  4 

az  A  2 


(5) 


where 


Ep  —  —  ( 1  +  f(~ )) 
•\P 


(6) 


is  the  grating  wavevector  of  the  p-th  chirped  grating. 

The  above  coupled  mode  equations  Eq.  (4)  and  Eq.  (5)  can  not  be  easily 
solved  analytically.  We  have  obtained  numerical  solutions  by  using  the  fifth- 
order  Runge-Kutta  method  with  adaptive  step  size.  Since  we  have  neglected  fast 
varying  phase  terms  in  deriving  the  coupled  mode  equations,  it  will  be  consis¬ 
tent  to  neglect  those  terms  with  very  large  phase  mismatch.  In  our  numerical 

:s 


calculations,  we  have  neglected  terms  with  phase  mismatch 


A  A-  =  I2A-  - 


>  0.5A-. 


-  ) 


With  this  cutoff  value  of  phase  mismatch,  the  program  runs  much  faster  than 
directly  integrating  Eqs.  (4)  and  (5). 


By  calculating  reflectivities  for  plane  waves  with  different  wavelength,  we 
c;m  obtain  the  reflection  spectrum.  Fig.  2  shows  the  reflection  spec  .nun  of  three 
chirped  gratings.  We  have  used  the  following  parameters:  n,i  =  1.5.  o  =  0. 


L  =  10/<m,  Ci  =  0,  02  —  f.  03  =  A]  =  0.5 /mi.  A_.  =  1.0/m/.  A.-j  =  1.4/m/, 
»i  =  0.02  4-  0.01(f)2  -  0.01(f),  n2  =  0.03  +  0.01(f)2  -  0.01(f).  // 1  =  0.01  + 
0.01(j)*  —  0.01(j).  e  i  -  — O.l(-jr),  (■>  —  O.l(j),  f;i  =  0.05(  j).  where  A  ;i  is  t  lit* 
center  wavelength  which  is  related  to  \p  by  AJ;  =  We  notice  that  by  chirping 
the  grating,  the  side  lopes  are  suppressed  and  the  band  width  is  increased.  In 
addition,  different  gratings  will  interfere  with  one  another  only  in  the  region 
where  they  overlap. 


3.  Filter  Design 

The  shape  of  the  reflection  spectrum  of  a  set  of  chirped  gratings  depends 
on  their  grating  amplitudes,  chirping  rates,  thickness  of  the  medium,  and  their 
relative  phases.  We  can  adjust  these  parameters  to  design  a  filter  with  desired 
reflection  spectrum.  By  increasing  the  grating  amplitude,  we  can  increase  the 
peak  reflectivity  and  the  band  width  of  the  reflection  spectrum.  However,  since 
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the  dynamic  range  of  the  gratings  is  usually  limited,  sometimes  we  cannot  reach 
the  required  reflectivity  by  increasing  the  amplitude  of  the  gratings.  Another 
way  to  increase  the  reflectivity  is  to  increase  the  thickness  of  the  medium.  A 
thicker  medium  can  give  a  higher  reflectivity  with  a  relatively  smaller  grating 
amplitude.  In  addition,  a  thicker  medium  will  also  result  in  a  smaller  band 
width.  However,  as  we  increase  the  thickness,  we  will  have  more  side  lopes. 
These  side  lopes  can  be  suppressed  by  tapering  and  chirping  the  grating.  On  the 
other  hand,  tapering  or  chirping  the  grating  will  also  increase  the  band  width  of 
the  spectrum.  Sometimes,  a  large  chirping  rate  will  produce  two  peaks  in  the 
reflection  spectrum.  When  we  use  more  than  one  grating,  the  relative  phases  of 
the  gratings  will  also  affect  the  spectrum  in  the  region  where  they  overlap.  In 
order  to  fit  a  required  reflection  spectrum,  we  need  to  consider  the  effect  due  to 
all  these  parameters. 

We  have  designed  several  filters  with  multiple  chirped  gratings.  Fig.  3  shows 
one  of  the  examples.  Fig.  3a  is  the  desired  reflection  spectrum.  Fig.  3b  shows 
the  reflection  spectrum  of  four  chirped  gratings.  V  e  notice  that  the  two  profiles 
match  quite  well.  The  parameters  used  to  generate  Fig.  31)  are  given  by  no  =  1.5. 
n  =  0,  L  =  10/m),  0\  =  0,  <f>2  =  7T,  —  0,  cq  =  0,  A i  =  0.315/m).  A_>  = 

0.33G5/)))).  A  )  =  0.3580/o)i,  A.,  =  0.39/m)..  »,  =  0.009  +  0.008(f)-  -  0.00S(f). 
))  >  =  0.0005  +  0.004(f)-  -  0.004(f),  n,  =  0.00S5  +  0.004(f)-  -  0.004(f),  m,  = 
0.018  +  0.007(f)2  -  0.007(f),  c,  =  -0.015(f),  f_>  =  0,  e,  =  0,  e,  =  0.015(f). 
This  example  demonstrates  the  reflection  spectrum  of  multiple  chirped  gratings 
with  relatively  large  separation  between  their  center  wavelengths. 

Another  example  is  shown  in  Fig.  4.  Fig.  4a  shows  the  two  desired  reflection 
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curves.  Fig.  4b  and  Fig.  4c  show  the  two  fitting  curves.  This  example  shows  the 


reflection  spectrum  of  multiple  chirped  gratings  with  their  center  wavelengths 
close  to  each  other.  By  adjusting  the  thickness  of  the  medium,  the  amplitudes 
of  the  gratings,  the  relative  phases,  center  wavelengths,  and  the  tapering  and 
chirping  factors,  we  can  fit  the  two  different  curves  in  Fig.  4a  by  using  two 
chirped  gratings.  Fig.  4b  and  Fig.  4c  show  the  fitting  curves  corresponding 
to  the  upper  and  lower  curve  in  Fig.  4a  respectively.  The  parameters  used  to 
generate  Fig.  4b  are:  r?o  =  1.5,  a  =  0,  L  —  10/mi.  6\  =  0.  =  ~,  Aj  =  0.9S5/m?. 

A)  =  1.05 urn,  n,  =  0.1  +  0.1(f)2  -  0.1(f),  n>  =  0.0S  +  0.04(f)2  -  0.04(f). 
<u  =  0.07(f),  e>  —  —0.05(f).  The  parameters  used  to  generate  Fig.  4c  are: 
no  =  1.5,  a  —  0,  L  =  22/mi,  <b\  =  0,  <?■>  =  0,  A]  =  0.95/mi.  A->  =  1.17/mi. 
n ,  =  0.03  +  0.01(f)2  -0.01(f),  n2  =  0.02  +  0.04(f)2  -0.04(f).  c,  =  0.01(f). 
e>  =  0.04(f).  We  notice  that  thicker  medium  and  smaller  grating  amplitude 
result  in  narrower  spectral  band  width. 


4.  Conclusion  and  Discussion 

We  have  analyzed  the  reflection  spectrum  of  multiple  chirped  gratings  using 
coupled-mode  theory.  Based  on  the  analysis,  we  have  designed  filters  with  desired 
reflection  spectrum  by  adjusting  the  parameters  of  multiple  chirped  gratings. 

In  the  above  discussions,  we  have  used  slowly  varying  envelope  approximation 
and  neglected  the  terms  with  large  phase  mismatch.  This  approximation  is  very 
practical,  if  we  need  only  to  know  the  general  shape  of  the  reflection  spectrum. 
When  we  design  filters  to  fit  desired  reflection  spectrum,  we  simply  adjusted  the 


parameters  according  to  the  past  experience.  A  more  general  approach  would  be 
to  vary  the  parameters  systematically  so  that  the  optimum  fitting  curve  can  be 
found. 
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Figure  Captions 

Fig.  1  Schematic  drawings  of  the  incident,  transmitted  and  reflected  beams. 

Fig. 2  Reflection  spectrum  of  three  chirped  gratings  where  r> q  =  1.5,  a  =  0. 
L  =  10 /mi,  0 i  =  0,  02  —  j,  03  =  A)  =  0.5/377?,  A_>  =  1.0/m?,  Aj  =  1.4/7777, 
771  =  0.02  +  0.01(f)2  -  0.01(f),  712  =  0.03  +  0.01(f)-’  -  0.01(f),  7?3  = 
0.01  +  0.01(f)2  -  0.01(f).  e,  -  -0.1(f),  e2  =  0.1(f).  and  f:!  =  0.05(f). 

Fig. 3  Reflection  spectrum  of  (a)  the  desired  profile,  and  (1?)  the  fitting  curve  with 
four  chirped  gratings:  ??o  =  1.5,  cv  =  0.  L  =  10/7???.  oi  =  0.  o>  —  rr,  0.3  =  0. 
04  =  0,  A 1  =  0.315/7777,  A)  =  0.3365/7777,  A 3  =  0.35S0/7777.  A.J  =  0.39/7777. 
771  =  0.009  +  0.008(f)2  -  0.008(f),  77-,  =  0.0065  +  0.004(f)2  -  0.004(f). 
773  =  0.0085  +  0.004(f)2  -  0.004(f),  774  -  0.01S  +  0.007(f)2  -  0.007(f), 
t\  =  -0.015(f),  e,  =  0,  e3  =  0,  e4  =  0.015(f). 

FigA  (a)  Reflection  spectrum  of  the  two  desired  profiles,  (b)  The  fitting  curve 
for  the  upper  curve  in  (a)  with  two  chirped  gratings:  ?7|j  =  1.5.  a  =  0,  L  = 
10/7777,  01  —  0,  02  =  7T,  A]  =  0.985/7777,  Xi  =  1.05/777),  77 J  =  0.1  +  0.1(f)2  — 
o.l(f),  772  =  0.08  +  0.04(f)2  -  0.04(f),  e,  =  0.07(f),  e,  =  -0.05(f). 
(c)  The  fitting  curve  for  the  lower  curve  in  (a)  with  two  chirped  gratings: 
770  =  1-5,  a  =  0,  L  =  22/77)7  ,  0i  =  0,  <*>>  =  0.  Aj  =  0.95/77??,  A 2  =  1.17/7???. 
77,  =  0.03  +  0.01(f  )2-0. 01(f),  ?72  =  0.02  +  0.04(f  )2-0. 04(f),  c,  =  0.01(f). 
e2  =  0.04(f). 
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Data  from  "case2.pl" 


cr 


X 


c*. 


X 

Lu  ^ 

o  55 
2  o 
o  o 

rz  co 

Q 
LU 

<  2 

^  Q 
<  — 


o 


o 

LU 

_J 

LU 

LU 

X 


LU 

Q 

LU 

X 

o 

H 

< 


o 

co 

< 

LU 

O 


o 

LU 


X 
LU 
h  LU 


Q 

LU 

X 

H 

X 

o 

LU 

X 

1“ 

o 

X 

LU 


Q.  ^ 

o.5 


X 

O 

l — 

o 

LU 
I — 
LU 
Q 


O 

O 


CO 

< 

0 


CO 

3 

Q 

X 

3 

o 

X 

0 

X 

o 

< 

m 

x 

i — 
X 
< 
LU 


< 

I — 
CO 


< 

0 

< 


CC  a) 


Atmospheric  model  used  to 
calculate  Earth  background 
radiance 


-&■  Mid-Latitude  Summer 


